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. Given N relativistic scalar free particles described by N mass-shell first 

(N 



class constraints in their 8N-dimensional phase space, their N-time descrip- 
tion is obtained by means of a series of canonical transformations to a quasi- 
Shanmugadhasan basis adapted to the constraints. Then the same system is 



. reformulated on spacelike hypersurfaces: the restriction to the family of try 



perplanes orthogonal to the total timelike momentum gives rise to a covariant 
intrinsic 1-time formulation called the "rest-frame instant form" of dynamics. 
The relation between the N- and 1-time descriptions, the mass spectrum of the 
system and the way how to introduce mutual interactions among the particles 
are studied. Then the 1-time description of the isolated system of N charged 
scalar particles plus the electromagnetic field is obtained. The use of Grass- 
mann variables to describe the charges together with the determination of 
the field and particle Dirac observables leads to a formulation without infinite 
self-energies and with mutual Coulomb interactions extracted from classical 
electromagnetic field theory. A comparison with the Feshbach-Villars Hamil- 
tonian formulation of the Klein-Gordon equation is made. Finally a 1-time 
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covariant formulation of relativistic statistical mechanics is found. 
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I. INTRODUCTION 



Relevant physical models are described by singular Lagrangians, so that their Hamilto- 
nian formulation is based on the Dirac theory of constraints M. The qualitative aspects of 
this theory are now understood , in particular the role of the Shanmugadhasan canonical 
transformation || in the determination of a canonical basis of Dirac's observables allowing 
the elimination of gauge degrees of freedom from the classical description of physical systems 
M. This programme was initiated by Dirac [H for the electromagnetic field with charged 
fermions. More recently, Dirac's observables for Yang-Mills theory with fermions (whose 
typical application is QCD) have been found || in suitable function spaces in which the 
Gribov ambiguity is absent. Also the ones for the Abelian Higgs model are known @ and 
those for the 577(2) x U(l) electroweak theory with fermions can be found || with the same 
method that works for the Abelian case. The main task along these lines will now be the 
search of Dirac's observables for tetrad gravity in the case of asymptotically flat 3-manifolds. 

The price for having only physical degrees of freedom is the nonlocal (and in general 
nonpolynomial) nature of the physical Hamiltonians and Lagrangians, as already known 
from Dirac's work on the electromagnetic field ||, for which the origin of the difficulties is 
the Coulomb self-interaction of the fermion fields. Two obstacles appear immediately: 1) the 
lack of manifest Lorentz covariance of the Hamiltonian formalism, which requires the choice 
of a 3+1 splitting of Minkowski spacetime; 2) the inapplicability of the standard methods 
of regularization and renormalization due to the nonlocality (and nonpolynomiality) of the 
interactions and the failure of the power counting rule. 

In Ref. ||, I reviewed these problems and I pointed out that the Lorentz covariance 
problem can be solved by reformulating the theory on spacelike hypersurfaces |lj] (see also 
Ref. U) and then reducing it to the hyperplanes orthogonal to the total momentum p^ (for 
the configurations in which it is timelike, p 2 > 0; see Ref. [10] for a general study of the 
embeddings of spacelike hypersurfaces in a given Riemannian spacetime). In this way, the 
breaking of Lorentz covariance is reduced to a minimum: only the three degrees of freedom 
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that describe the canonical centre of mass 3-position of the overall isolated system are not 
covariant. In Ref. ||, it is also stressed that only one tool has until now emerged for attacking 
the difficult task of quantizing a nonlocal (and generically nonpolynomial) theory. Namely, 
the problem of the center-of-mass extended relativistic systems in irreducible representations 
of the Poincare group with P 2 > 0, W 2 = —P 2 S 7^ (they are dense in the set of all allowed 
field configurations) identifies a finite worldtube of noncovariance of the canonical center of 



mass, whose radius p = \ / —W 2 /P 2 = \S\fvP 2 represents a classical intrinsic unit of 
length that can be used as a ultraviolet cutoff at the quantum level in the spirit of Dirac and 
Yukawa. As noted in Ref. [6], the distances corresponding to the interior of the worldtube are 
connected with problems coming from both quantum theory and general relativity: 1) pair 
production occurs if an attempt is made to localize particles at these distances; 2) relativistic 
extended bodies with a material radius less than p cannot have the classical energy density 
positive definite everywhere in every reference frame, and the peripheral rotation velocity 
may be higher than the velocity of light. Therefore, the worldtube is the flat remnant of the 
energy conditions of general relativity; in this theory, the radius p is defined in terms of the 
asymptotic Poincare group that exists in the case of asymptotically flat 3-manifolds. 

In this way, one has perceives the possibility of formulating a "rest-frame field theory" 
on these special hyperplanes ||, in which the ultraviolet cutoff just identified could be used 
in a constructive way. The asymptotic states of this theory should be connected with the 
description of free particles on these hyperplanes. But this implies that the phase space 
description of the particles involves only 6 degrees of freedom per particle, that there is an 



1-time description ( "rest-frame instant form" of the dynamics in the language of Ref. [11 
the time parameter being the rest frame time, which describes the foliation of Minkowski 
spacetime with these hyperplanes), and that there is no room left for relative times and 
energies, which are precisely the variables connected with the spurious solutions of the 
Bethe-Salpeter equation [[12] for relativistic bound states in ordinary quantum field theory. 
Actually, the standard Fock space describes asymptotic states defined as tensor products 
of free one-particle states without any restriction on their mutual temporal ordering (an 



in-particle may be in the absolute future of another in-particle) : the classical background is 
an N-time description of N free particles, each one of which has 8 degrees of freedom and 
an associated mass-shell first class constraint (connected with the inverse of the standard 
propagators) . 

Therefore, in Ref. || there was a review of the problems of relativistic particle mechanics, 
especially the two main ones: i) the No-Interaction Theorem (see Ref. J15 for a 



review), which turns out to be also present in Galilean mechanics |T(| and to be connected 



with the use of an N-time description; ii) the many definitions of relativistic center-of-mass 



position (see Ref. JT/J for a review), since the Lorentz signature of Minkowski spacetime 
precludes the existence of an object with all the properties of the nonrelativistic center-of- 
mass. 

Behind all these descriptions there is the theory of the irreducible representations of the 
canonical realizations of the Poincare group |18[ in given phase spaces. 

In this paper, I will consider N free relativistic scalar particles, and I will study their 
description both in the N-time approach (following jointly Ref. [TP]] and the complete analysis 
of the N=2 case in Ref. J20|) and in the 1-time approach, which was developed in Ref. [f2T] 



for completely different reasons. While the N-time theory identifies a set of variables that 
makes it possible to disentangle the relative times and energies from the relevant variables, 
the 1-time theory allows the identification of the branches of the mass spectrum of the N- 
body isolated system. Then the two descriptions are compared, and it is found that the final 
physically relevant variables are the same in the two approaches and that they also seem to 
be the best relativistic kinematical variables (carrying the knowledge of the geometry of the 
timelike Poincare orbits) for the study of interactions. 

In Section II, there is a study of the canonical transformations needed to find a quasi- 
Shanmugadhasan canonical basis adapted to the N-l combinations of the original first class 
constraints that define the vanishing of the relevant N-l relative energies. Due to the nonlin- 
earity of the final canonical transformation, it is not known how to find the mass spectrum 
and the inverse of the transformation. 



In Section III, the system is reformulated on spacelike hypersurfaces, on which each 
particle is now identified only by three coordinates. After the identification of the first class 
constraints, the reduction to the special intrinsic family of spacelike hyperplanes orthogonal 
to the total four- momentum, when it is timelike, is studied. The 1-time covariant rest-frame 
instant form of dynamics is identified. 

In Section IV, the N- and 1- time descriptions are compared, and the mass spectrum is 
found. 

In Section V, there is a study of the mass spectrum in the N-time theory. Then action- 
at-a-distance interactions are introduced among the particles in the 1-time theory with a 
comment on their separability. 

In Section VI, the isolated system of N charged scalar particles plus the electromagnetic 
field is formulated on spacelike hypersurfaces and then reduced to the hyperplanes orthog- 
onal to the total timelike four-momentum. The charges of the particles are described in a 
pseudoclassical way by means of Grassmann variables. The Dirac observables of the parti- 
cles and of the electromagnetic field with respect to electromagnetic gauge transformations 
are found. The final resulting four first class constraints contain the interparticle Coulomb 
potential (extracted covariantly from the classical electromagnetic field theory) but not the 
classical electromagnetic self-energy due to the vanishing of the square of the Grassmann 
charges: the underlying hypothesis of charge quantization generates a regularization of the 
pseudoclassical description. 

In Section VII, the 1-time theory of one charged scalar particle plus the electromagnetic 
field is compared with the Feshbach-Villars Hamiltonian formulation of the Klein-Gordon 
equation in an external electromagnetic field. 

InSection VIII, after a review of the main problems of covariant relativistic statistical 
mechanics, its reformulation in the 1-time theory is given. 

The Conclusions IX contain some final remarks and the outline of future research. 

In Appendix A, there is a review of the properties of the standard Wigner boost for 
timelike Poincare orbits. In Appendix B, the formulas for the case N=2 are reviewed. In 
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Appendix C, some notation in the description with spacelike hypersurfaces is given. 
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II. THE N-TIME THEORY 



A system of N free scalar relativistic particles is usually described in phase space by 8N 
variables xf,pf, i=l,..,N, {xf,p^} = —5^1]^ \r\^ v = (1,-1,-1,-1)] restricted by N first 
class constraints [to start with, we assume arbitrary masses m; for the particles] 

l=p 2_ m 2^ O; {0 i ,0 i } = O. (1) 

The evolution in a scalar parameter r is described by the Dirac Hamiltonian 

N 

H d = Y,Ht)<I>u (2) 
i=i 

where X%{r) are Dirac's multipliers. This description is obtained by starting from the La- 
grangian 

N 



C = -Y^m^xfir), (3) 
i=i 

whose associated canonical momenta are 

dC x^ 

Pi = -77— = m i~h- (4) 
aa V J x~ 

Therefore, the configuration position variables, describing the worldlines of the particles, 
are q% (r) = xf(r). An alternative N-time Hamiltonian description is the multi-time one [|J, 
in which the N scalar time parameters t% are defined by dri = Ai(r)dr. In this description one 
has gf(rj) = x% (rj), {xf(rj),pj (r,)} = —bijff v ', and the first class constraints 4>i ~ are the 
Hamiltonians for the evolution in the Tj's [ a ^g^.' Pfc ^ ={^4, (pi} are the many-time Hamilton 
equations and {4>i,(f)j} = are their integrability conditions; = means evaluated on the 
solutions of the equations of motion]. The Lagrangian description derives from the action 

s = e£i / dr i £i witl1 A = ~ m i J^K^i). 



Let us remark that in presence of interactions, the No-Interaction-Theorem [13,14] im- 
plies gf(Tj) 7^ (n, .., rjv). As shown in Ref. |I| (see also Ref. ||), this result is also 



present at the nonrelativistic level: it is not connected with the Lorentz signature but with 
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the many-time description. If the particles are charged, the minimal coupling to an external 
electromagnetic field requires the canonical covariant coordinates xf s and not the configu- 
ration 's, which describe the particle worldlines: therefore the electric charges cannot be 
localized on the worldlines in the interacting case. 

In Ref. jBJ, a study of the kinematics of N relativistic free scalar particles was started 



with a double aim. Firstly, the spurious solutions of the Bethe-Salpeter equation [12| stim- 
ulated a search for nonlinear canonical transformation replacing N-l linear combinations 
of the first class constraints fa ~ with suitable "relative energy variables". Secondly, it 
was investigated how to introduce act ion- at- a- distance interactions while preserving the first 
class nature of the constraints (only nonseparable interactions were found). 

Subsequently, for the case N=2 with p 2 > [p^ = Y^iLiPi is the total momentum] 
suitable canonical variables adapted to the timelike orbits of the Poincare group were found. 
Then, in Ref. [^j, a study of the N=3 case was initiated, with the aim of understanding 
how to describe separable action-at-a-distance interactions. 

Collecting all these results and simplifying the notation, we start with the canonical 
tranformation from the basis xf,pf, i=l,..,N, to a first set of linear center-of-mass and 
relative variables [a=l,..,N-l] 



X; 



(5) 



defined by 



= ViVEill laiXl {R%, Oft = - W 

where 7, 7 a is a basis of orthogonal unit vectors satisfying 

N N 



(6) 
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N N 

i=l i=l 
N-l 

lilj + E ^".iHaj = Sij. (7) 
a=l 

From now on, we shall choose 7, = l/y/~N, i=l,..,N, so that we have 

Y n N-i y 

^ = E ^ E 7ai7« = - Tr- ( 8 ) 

1=1 0=1 

Here x^ is a canonical covariant center-of-mass coordinate, which however does not have 
free motion, since = {x^,H D } = — -j| J2f=i K(j)p% / \f$ 7^ A(r)p M / 'ri^fp 1 = A(r)w M (p) 
[A(t) = ?7 = signp ]; it moves with a "classical zitterbewegung" that depends on the 
choice of the arbitrary Dirac multipliers \(r) (it is a gauge effect). 

The inverse canonical transformation is 

1 N-l 

< = *r + -7= E %iK 

1 _N-1 

= -I? + VN %iQl (9) 
iV a=l 

The conserved Poincare generators are 

N 



1=1 

V 

^ = E - = ^ + ^ 



i=l 



IT = x^tf - x v p» 

JV-l 

^=E(W-W). (10) 

a=l 

The isolated system of N particles is assumed to belong to an irreducible timelike repre- 
sentation of the Poincare group with Casimirs p 2 > and W 2 = —p 2 S [W* 1 = \^ VOL ^p v S a ^ 
is the Pauli-Lubanski fourvector (e 012 3 = 1, e u ' fc = e^); S is the rest-frame Thomas spin 
(see later on)]. 

Note that U J,V and are not separately constants of the motion due to the "classical 
zitterbewegung of x M " (which is an analogue of the zitterbewegung associated with the Dirac 
position operator for the electron). 
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Extending the results of Ref. |20| for N=2 to arbitrary N (see Appendix A for the 
notation), one defines the following further canonical transformation 





x" 




x fJ, pil 






T Ra 


R a Qa 




Pa 7? a 











(11) 



{x»,jT} = -rT, {T Ra , e Rb } = -S ab , {pi 4} = 8 ab S ij , 
with (one has p^ = 



gup 



vVi^ip + vVp 2 ) 



p'2 



(12) 



f- = p» 



T Ra = e;(u( P ))RZ 
e Ra = e°Ju(p))QZ 



P ■ Ra 



VVP 
P-Qa 

r}\fp 1 



P\ = e'(«(p)X = R r a 



P 



7]\JP 

■J" 



P- Ra 



P° + T\\TP^ 

"i]\Jp p + r]\P 



(13) 



where Eqs. ([A10|) and ( |A11|) have been used to evaluate x^. 

As shown in Ref. [^(J and in Appendix A, the canonical transformation of Eqs.(|l"3"D 
is defined by boosting at rest the relative variables R%, with the standard Wigner 
boost L A n{p,p) = ej(u(p)) of Eqs.(|A2|), Q for timelike orbits p 2 > 0. In Eq.flll]), 
T Ra = e°{u{p))R% and e Ra = ^(u(p))Q^ are Poincare-scalar relative times and energies, re- 
spectively; p r a = e r ^{u{p))R^ and n r a = e r (u(p))Q% are, respectively, relative three-coordinates 
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and three-momenta, which transform as Wigner spin-1 3- vectors. It turns out that the new 
canonical center-of-mass coordinate x M is not a fourvector under Lorentz boosts (it has only 
Euclidean covariance), so that it cannot define a geometrical center-of-mass worldline, but 
only a frame-dependent line (pseudo- worldline). 
The Lorentz generators become 



2 p dp^ dp u 



[p p (S^p u - S Pu p») + r]Jp~ 2 (S otl p u - 



V^/p 1 



[{P - VyP ) S + 



1 



S %1 = S tJ + —=(S ol p> - S OJ p % ) - 



P° + Tjy/p 1 

P k (s ki p> - 5 fc V) 



rj^/p 2 rj^/pFij) + rj^/p 1 ) 

If we introduce the rest-frame spin tensor (the last line gives the Thomas spin) 



(14) 



S AB = 4{u{p))e»{u{p))S^ 



N-l 



a=l 

s rs = E(pX-pX) 



a=l 



S r = le rst S st =J2S r a , S a = p a xn a 

1 a=l 



we find the following form of the Poincare generators 



(15) 



r = xy - x ] p i + s ir 5 js s rs 

5 ir S rs p s 



J m = x°p l - x l p° - 



p° + f]\fp z 



(16) 
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and the following form of x^ 

1 



x» = x" - 



fJ, ( qoA 



S Ar p r 



+ p^ + l-q^fpp-rf 

po _|_ rj^/pZ Tj^/jFij) + rj^fp 1 ) 



S 5r p r 



Y N-l 

X° = X° -J2 ( T RaP- K a ~ e Ra p- p a ) 

P a=l 



N-l 



X % = x l 



VVP 2 a =l 



PaP ■ ~ KP ■ Pa , P l {T Ra p ■ lt a - t Ra p ■ p a ) 



P° + T}\fp l 



+ 



r]^fp 1 {p° + rj^fp 2 ) 



Therefore, we have 



5 ir S rs p s 



s m = 



P° + Tj^/p 1 



= S AB - ^M^MM - e B MP))^ KiUip)h ° 



iAB 



dp/i 



e a D (u(p))r) EC S 



CD 



The inverse canonical transformation is 
1 



JV-l 

X° = X° + Yl ( T RaP- - £RaP- Pa) 
P a=l 

N-l 



X = X + 



1 



+ 



VVP 2 a =l 
1 



Y (T Ra 7i a - e Ra p a ) + 



N-l 



7=t7 - , Y [P • P<^a ~ P- T^aPa H T^{T R aP '• - ^RaP 

VP{P +VVP)a=l VVP 



r: 



i 



-AT^ + P-Pa) 



tia - Pa^ r^Uft. + — 7=rJ 

r]^p z p° + r]^/p z 



Q°a 



1 



V^/p 1 



{£RaP° + P- 7T a ) 



A P , P-^a v 

Q a = n a -\ i=^{e R a + — 1=^) 

ri^/p 1 p° + 77VP 
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and the original variables xf , pf have the following form in this new canonical basis: 



x° = x° + 



x n-1 x l 

V] [-^^lai{T R aP° + P- Pa) H ^(^W ' ~ ^RaP Pa)} 

— ' ' A r/VP 



?7vV ^ ViV 

JV-1 



1 VfT - -> \ I V- Pa^a-V^aPa . 



I — T / j L rta"a *-naraj < „ , / — ?r 

WP a=l P + WP 

P TRap- TX a - e Ra p- p a 1 A_J 



+ 



P° ViV ^ 



V AT ~[ WP P + WP 



P-P a 



wp 7 p° + vVp 1 

N-l 



N WT t^i 



Vi = j^ + ^NY, 7a 4 [vr a + -^j^ (e Ra + — 



N-l 



P 



P ■ 7T a 



0=1 



W^* 7 P° + f]\fp 1 



)]■ 



(20) 



As shown in Ref. ||20|| , it is then convenient to perform the following canonical transfor- 
mation 



defined by 









T 








z k 






T Ra 






e 


Pa TTa 






Pa 7?a 


T Ra 








£ R a 



e = r\\lp* 

p-x p-x 
-t = — f=f = — f=s = m(p) • 2; 



WP 2 WP 2 ^ 



fc = u(p) 



p 



WP^' 



k° = u °(p) = yi + k 2 



Z — ri\frp(x — —x°) 



(21) 



= {T,e} = -1. 



(22) 



Its inverse is 
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e e 



p° = e^l + k 2 

p = ek, (23) 



and the form of the Poincare generators is now 



jn = z i k j _ z j k i + 6 ir 5 js S rs 

sir crs n s 

j oi = - z *k° + v J io = -J m . (24) 



e(l + \Jl + k 2 ) 

Here T is the Lorentz-scalar rest-frame time and e is the invariant mass of the isolated 



system of N particles. Actually, Eqs.(21) define as many canonical transformations as there 
are disjoint branches of the mass spectrum (half with rj = signp > and half with rj < 0). 
All the noncovariance of x M has been shifted to the 3- vector z [it has the Euclidean covariance 
(see Appendix A) of the kinematical stability subgroup of the Poincare group of timelike 
orbits]. The noncovariance of and i? is a consequence of, also in the free case, the Lorentz 
signature of Minkowski spacetime and of the geometry of timelike Poincare orbits, when 
one uses the associated instant form of dynamics according to Dirac ||11|| . In contrast to 



what is often said, it is not due to the relativistic version of the No-Interaction-Theorem: 
until it was recognized that this theorem is independent of the Lorentz signature jy|, it 



was difficult to separate this phenomenon from the complications of relativistic kinematics. 
In the canonical subspace z, k (treated independently from the canonical pair T, e, whose 
presence is due to the covariant description with first class constraints) one has z = rjy/p 1 ^, 
where the canonical noncovariant three-coordinate k = —K/p° + 5 x p/p°(p° + rjy/p 1 ) [from 
Eq.(p4|) with K % = J m ] can be shown to depend only on the Poincare generators (see 
the theory of the canonical realizations of the Poincare group |18|] as a special case of the 



canonical realizations of Lie groups [23]), to be the classical basis of the Newton- Wigner 
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three-position operator [24|, and to have free motion without classical zitterbewegung |25 



(the problem of classical zitterbewegung will be clarified in Ref. p5|). While these results 
were obtained in the standard instant form of dynamics x° = const., now they are recovered 
(in terms of z) in a covariant version of the instant form, the "rest-frame instant form 
T = const. 11 . The vector z describes the Jacobi data for the center-of-mass three- variable 
in this rest-frame instant form and J^ u in Eqs. (|24]) is split in two terms that are separately 
constants of the motion; in this respect, x M has all the properties, except covariance, of the 
Foldy-Wouthuysen mean position for the electron |27]] (see also Ref. [ 25| , where a covariant 



canonical zitterbewegung-free x^ is shown to exist for certain pole-dipole systems, though 
no statement is made about extended systems). As shown in Ref. |p5|,|£|,|9l , there are other 
two notions of center of mass that can be build solely in terms of the Poincare generators, 
coincide with k in the rest frame and have the same three-velocity of k. In the instant form 
x° = const., they are the Fokker center of inertia |28| Y and the Moeller center of energy |29| 



R; in the rest-frame instant form T = const, they become ny/jfiY and rjy/j^R, respectively. 
The vector Y is covariant but not canonical ({Y\ Y^} ^ 0), while R is neither covariant nor 
canonical. Since in the rest frame H = Y = R, in every frame the Fokker center of inertia 
describes by construction the worldline obtained by applying to the rest-frame k the Lorentz 
transformation appropriate for going from the rest frame to the given frame. In contrast, 
R is defined by replacing the masses n%i in the Newtonian center of mass by the energies 
p°. Both k and R define in every frame a pseudo- worldline; if we draw all these pseudo- 
worldlines in a given reference frame, we get a worldtube around the worldline of the 



Fokker center of inertia, whose scalar intrinsic radius p [29] is determined by the Poincare 
Casimirs for timelike orbits with spin 



V VP 

See Refs. |||9|] for the remarkable properties of the worldtube and for the proposal to use p 
as a ultraviolet cutoff. 

It should be remarked that the definitions of k,, Y, R are naturally given in phase space. 
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Even in the free case (gf = xf), their definitions in configuration space in terms of the gf s 
are very involved, being velocity-dependent. This is the reason why the Newtonian definition 
x = Y,iLi m iQi/ J2iLi m i m configuration space does not extend trivially to the relativistic 
case. 

However the canonical variables of Eqs.(|2"T|), ( P2"D are not the final ones adapted to the 
first class constraints of Eq.(fj]). This can be seen by replacing those constraints with the 
following linear combinations 

N N N-l 

X = NY / ^=P 2 - NY\m* lultiQc ■ Qd = 

i=l i=l b,c=l 

N N-l N N-l 

= P 2 -N(Z m^-Nj2Ql)=e 2 -N[j: m\ - N £ {e\ b - vr 2 )] « 

i=l 6=1 i=l 6=1 

. r, r 2 AT N ^\ A „ 

Xa = — 2^ JaiPi =P'Qa — 2^ 7aiK ~ ^ 2^ IbilciQc ' Qc\ = 

A i=l A i=l 6,c=l 



ifNtf* N N - 1 
A P i=l b,c=l 

\/~N N N ~ 1 
= ee Ra — laA m2 i ~ N %aci(e Rb e Rc - 7i b ■ n c )] = ee Ra « 

1 i=l b,c=l 

1 2 N - 1 1 2e N ~ 1 

& = JpX +^1Z ^iXa = J^X + E ^Ra « 0, (26) 

where e^ a will be explicitly defined in Eqs. (|3l|) later on. 

We see that the N-l constraints Xa ~ do not determine the e^ a directly, but they form 
a system of N-l quadratic algebraic equations for them. With the constraint x ~ we 
obtain a system of N equations for e and the e^ a 's. By substitution we could arrive at a 
single (e-dependent) equation for each e Ra of order 2 7V_1 ; therefore, in general there will be 
2 N ~ l sets of different solutions for the e# a 's, a=l,..,N-l. Each one of these sets of solutions, 
substituted into x ~ 0, should give the equation for the spectrum of the invariant mass e 
of the isolated system of N particles in terms of m 2 and 7T& ■ tt c . Since each free particle has 
two disjoint branches of its mass spectrum (i]i\fp^ ~ the whole system will have 2 N 

disjoint branches for e (2 N ~ 1 for p 2 = e 2 ); 2 N is a topological number, namely, the dimension 
of the zeroth homoty group of the mass-spectrum hypersurface in phase space. Even in the 

17 



free case, these branches are not known because of the difficulty of solving the nonlinear 
equations for the e^ a 's. 

By defining H D = X(r)x + J2a=i A a (r)xa, one sees that the splittings J"" = L^ v + = 
L^ v + S^ v do not define separate constants of the motion I/ a/ , S^ u or L^ u , S^ v ; only Eqs.(^) 
split J^ v into a pair of constants of the motion. 

In any case, Eqs.(^) suggest the following nonlinear (point in the momenta) canonical 
transformation |1| of the canonical variables of Eqs.@ 



X^ £> M 



(27) 



N 



N-l 



x 1 - 



& + —(rT - -^f-) EK - N E luaMu ■ Q v ) x 

P i=l u,v=l 



2p 2 

N-l 



x i-nrf E 1™ R ™ - 
ViV a=1 

N-l N 
— ^ P ' RbQ 'c^E T v A'v^^-i /"I 

6,c=i i=i P + N\/N J2 d=l ldj p ■ Q d 



{N5i, - 1)7, 



x (ibj + 



u=1 luhP-Qu 

^ h=1 pi+NVNj:^ 1 % hP -Q v 



)} 



N 



N 



N-l 



X h 



i=l 



u,v=l 



N-l 



J2lai(e?(u(p)) P r a -u»(p)T Ra ) 



E T Rb ^{u{p))< - ^(p)e Ra ) E - ^ ~^? C , 

=1 e + NVN Ed=i 7*-£Ai 



,,c=l 



x (tw + 



EN-1 _ 



)] 



P^ =P^ 



N-l 



N 



R£ = R£- nVN E P • #&Q?E 



l,C=l 



x 
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EN-1 « _ 
u=1 lukP-Qu 

= R»-nVN V T^Q(fV 7=^-1 

EN-l „ 

rM, N N-l 



Q'a = Qa 1 ~ ^^J2%i(mt - N " £ %laQ b ■ Qc 

A p 



i=l 6,c=l 



= Qa 9 5" J2 7ai [m? - iV ^ luld(£Rb£Rc ~ n • 7T C )] • (28) 

/ P i=l 6,c=l 

We now have [it can be checked that eqs.fllOD are reproduced by inserting i?^, Q% 
from Eqs.(ES|)] 



iV-l 



a=l 



The subsequent canonical transformation 







pfJ. 








TRa 






Pa K a 











defines the following analogues of the variables of Eqs. ([l3|) 
X =X + 2 6 "^ P ^ VAB ~ dp 





= F - 
7 




= jf 




_P- Ra 




r}\fp 1 



1 



p 



N-l 



(29) 



(30) 
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EJV-1 . 

^ lailci ^ Z " fc=1 e+NVN^-i WRv v, 

. _ P-Qa _ 1 
C_Ra / — rr Xa 

1 v^v N N ~ 1 
= ~\ ee Ra 7T~ ^Zlaiijnl - N 7w7«[ej»ejte - ^6 • ^c])] ~ 

6 2 i=l b,c=l 

_ N-l 

X = Pa~ NVN ]T T m 7T c X 
6,c=l 



v-^JV-l . 

y y 7a.7« ^^WlvE^ 1 ^, 
i=1 e + ivy iV Ed=i 7d»ejM ^ fe£fl " 

7?a = 7Ta- (31) 



We now have 



jv = L +S 



L = X p — X £T 



S =S ^- - e ^u(p)) VAB 



de»(u(p)) ^ de»(u(p)) 



dp, 



■P ~ 



dp u 



p»)S pa 



1 



f]\fp 1 {p° + Tj^/p 1 ) 



[pp{S^p u - S Pu p») + r]Jp~\S otl p u - S ov p»)} 



~ oi 1 r- - ■ v k (S ko v i - S ki v° 

S = -^[(p°-r l Jf)S m + P 1 P ' 



p° + ri^fp 1 
p k (S ki p j - S*V) 



s 3 = S* + -±= (S°Y - S°ip l ) - _ . 

?7VP r]^p z {p° + r]^p z 



S = e*{u{p))e B y {u{p))S^ 



a=l 

= zVx - pX) 

a=l 
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5 = -e rst S = $a, S a = p a x n a 

/ 0=1 



-u 1 -oA 5 r> r 



H -^b p 



ri^pp 1 p° + i]\fp 1 •\]\fp 1 {j) + •\]\fp 1 

Note that if we put Tr u = 0, a=l,..,N-l, then we get 



(32) 



■ Ra 



0. 



Pa = Pa, 



s =s l 



N-l 



r/^S p r 



a=l 



X ~ X 



p-5=p-x = p-x. 
The final canonical transformation, the analogue of Eq. (|2"TD , is 









f 


X 






z k 






TR a 






e 


Pa ^a 


(-Ra 




Pa Ka 


TRa 








(Ra 



with 



T 



p ■ x p ■ x 



ri^fp 1 ri^fp 1 

2 



e = rj\ip 

z = rjdp A [x x J 



p" 



k 



P 



rj^/p 1 



(33) 



(34) 



(35) 



The canonical basis of Eqs.(|3~4]) is a quasi- Shanmugadhasan basis because the N-l relative 
energy variables vanish if the constraints \ a ~ are used; therefore, it seems suitable for 
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the description of interactions, which yield relativistic bound states without relative energy 
excitations at the quantum level. Only the mass-spectrum constraint is not contained in the 
basis, because it describes 2^ disjoint branches 

e- X A (mi,n a -n b ) « 0, A = 1,..,2 N . (36) 

Therefore, in the free case one could define 2 N Shanmugadhasan canonical transforma- 
tions containing also the variables e — A^; this is also possible for all Liouville integrable 
interactions. Moreover, Eqs. (|34D are adapted to the Poincare Casimir p 2 = e 2 , single out a 
rest-frame time T as the natural time parameter, and show that the relative times T Ra are 
gauge variables: namely, it is a freedom of the observer whether to describe the N particles 
with a 1-time theory (for instance by adding the gauge-fixings Tr„ = 0) or with an N-time 
theory [or even with an F-time one with 1 < F < N]. 



The only problem with the canonical transformation ([H]) is that to get its explicit in- 
version one has to solve the system of equations 



2 

^ 6,c=l i=l 



N N-l 



laA m2 i + N J2 Ibilcrfb ■ 7T C ] = 0, (37) 

1 i=l 6,c=l 

in order to obtain the old variables e# a in terms of the new e# a and the kinetic terms ttj, ■ n c . 
For e^ a , this is the same set of equations discussed previously, whose solution is needed also 
to determine the branches of the mass spectrum. So far, the inversion is known only for 



N=2 poyi9| , as can be checked in Appendix B. Actually, it turns out that one should solve 
Eqs.([37|) together with the first of Eqs.(f2~6|), 

X = e 2 - iV[£f =1 mf - NEt'^ - vf 2 )] » 0, 
which determines the mass specturm, as a system of N equations in the N variables e, e^ a , 
a=l,..,N. 

Equations (|37|) are still untractable even if we restrict ourselves to a rest-frame 1-time 
theory by adding the gauge-fixings T Ra w 0, which imply T Ra « and all the results of 
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Eqs.([33|), so that one can eliminate the N-l pairs of second class constraints T^ a « 0, 
tRa ~ by going to Dirac brackets. 

We shall see in the next Section that there is an approach to the 1-time theory that 
allows us to find the 2 N branches of the mass spectrum and the 2 N ~ 1 solutions of Eqs.(p7|). 
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III. THE 1-TIME THEORY 



Since in the case of a charged scalar particle interacting with the electromagnetic field it 
is not clear how, in the absence of a covariant notion of equal times, to evaluate the Poisson 
bracket of the particle mass-spectrum constraint (p(r) — eA(x(r))) 2 — m 2 m and of the 
field primary constraint tc°(z, z°) ^ 0, this system was reformulated in Ref. |[T| on spacelike 



hypersurfaces following Dirac's approach to parametrized field theory PJ. To describe the 
intersection of a particle worldline with S(r), only three coordinates a = i/(r) are needed 
and not four, since the time variable is the parameter r labelling the hypersurfaces of the 
family. But this implies a covariant solution of the constraint p 2 — m 2 »s and a choice of 
the sign rj = signp°. For N free particles, one obtains a 1-time description for every choice of 
rji = signp° with particle coordinates fji(r) and with (r) = ^ M (r, r/j(r)) describing the i-th 
particle worldline [so that if (r) = z£(t, rji{r))+Zf{r, f}i( T ))Vi ( T )]'i moreover, signx°(r) = rji. 
See Appendix C for the notations. 

As shown in Ref. [21|, the system is described by the action 

S = J drd 3 aC{r,a) = J drL{r) 

N I : 

C(t, a) = ~^2S 3 (a - r/i(r))^m iV / 'g rT {r,a) + 2g Tf (r, a% r (r) + g fS (r, ^[(t^t) 

i=l 

N 



L ( T ) = -J2Vi m i\/9Tr(r,f!i(T)) +2g T f(T,ff i (T))r]l(T) + gfs(r,rfi( T ))Vi( T )Vi(^), (38) 

i=l 

where the configuration variables are z m (t, a) and ^(r), i=l,..,N. The action is invariant 
under separate r- and a-reparametrizations. 
The canonical momenta are 

. dC{r,a) o 

Pn\T, a) = - = 2^ 5 {a- ^(r^m; 

UZ T [T, a) i=1 



^g TT (r, a) + 2g Tf (r, ^f(r) + g f g(r, a)r)[(r)?)|(r) 
[(pjni,+ (PuZ^)Y s z s ,}(r,a) 
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= Vi m i I = ; (39) 

y9rr(r, ffi(r)) + 2g Tf (r, ^(r))?)[(r) + g„(r, Vi(^))Vi (r)v! (r) 

and the following Poisson brackets are assumed 

{z^T,a),p u (r,a'} = -^5 3 (a-a') 

K(r),/c}(r)} = ^*. (40) 



From Eqs.(39) one obtains the following four primary constraints by using Eqs. 

sr 



N 



W m (t,ct) = p M (r, a) - i M (r, 5) ^ <5 3 (ct - ^(r^ymf - j fS (r, o^if ( r )M r ) - 

i=l 
iV 

- ^ M (r, a)7 r ^(r, a) £ 5 3 (a - ^(t))k« « 0, (41) 

i=l 

which satisfy 

{^(r,a),^(r, ( r')} = 0. (42) 

Since the canonical Hamiltonian vanishes, one has the Dirac Hamiltonian [A M (r, tr) are 
Dirac's multipliers] 

#d = / dVA^^O-.o), (43) 

and one finds that {Ji^r, a), Hp} = 0. Therefore, there are only the four first class con- 
straints of Eq.fl4~l|). The constraints 7Y M (r, <r) w describe the arbitrariness of the foliation: 
physical results do not depend on its choice. 

The conserved Poincare generators are (the suffix "s" denotes the hypersurface S(r)) 

p» = Jd 3 ap»(r,a) 

JjT = I d 3 a[^(r, t)p u (r, a) - z»(r, a)ff(r, a)}, (44) 



and one has 
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{z!>(T,d),fi} = -rfi (45) 

r N i 

/ dVW M (r, B) = - P{t, ffiir^^ml - t«(t, ^(r))^(r)^(r) - 

J i=i 

- E^^W)W.^W)««(r) w 0. (46) 



i=l 



Let us now restrict ourselves to spacelike hyperplanes Eff(r) by imposing the gauge- 
fixings 



"(r, a) = z"(r, a) - a^(r) - #(t)</ « 



{^(r,a),^(T,a')} = - V Z8 3 (B-B'), (47) 

where fe^(r), f = 1, 2, 3, are three orthonormal vectors such that the constant (future point- 
ing) normal to the hyperplane is 

p(r, B)^F = b» = e^6?(r)&f (t)^(t). (48) 

Therefore, we get 

z£(t, B) « #(t) 
2f(T,3)«if(T) + ijf(r)a' 

<7„(t, a) « -<J„, 7 «(r, a) « -5", 7 (r, a) « 1. (49) 
By introducing the Dirac brackets for the resulting second class constraints 
{A, B}* = {A, B}-J d 3 a[{A, C(r, B)}{H,(r, B),B}- {A, H,(r, B)}{C(r, a), B}], (50) 
one finds by using Eq.(pE5|) 

W(T)MT)}* = -rr. (51) 

The ten degrees of freedom describing the hyperplane are x%(t) with conjugate momen- 
tum and six variables <p\ (r), A = 1, .., 6, which parametrize the orthonormal tetrad b^(r), 
with their conjugate momenta T\{r). 
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The preservation in time of the gauge- fixings (^(t, <t) ~ implies 

^C(r, B) = {C(r, ff), H D ) = -A"(t, a) - ^(r) - &£(r)</ « 0, (52) 

so that one has [by using = and 6^(r)6^(r) = — &£(r)6£(r)] 

A"(r, a) «A"(r) + n(r)^(rK 
A"(r) = -±?(r) 

A^(r) = -A^(r) = ^(t)^(t) - ^(r)^(r)]. (53) 
Therefore the Dirac Hamiltonian becomes 

H D = A"(t)W„(t) - h^(r)H^(r), (54) 
and only the following ten first class constraints are left 

N N 



W{r) = f d 3 aW(r, a) = p» - Z" £ r^mf + /f (r) + &J?(r) £ ^(r) « 0, 
J i=i i=i 



N 



= ST to - mr)K - bt(T)K)Y,rt(T)rii\/rn2 + «f (r) - 

i=l 

AT 

- P£(t)6&(t) - 6£(t)6S?(t)] £<(tK(t) « °- ( 55 ) 

i=l 

Here 5^ is the spin part of the Lorentz generators 

Jg = XgP s — x s Pg + Sg 

S-f = #( r ) | Saa f p v {r, a) - b v f (r) J ^/(r, a). (56) 



As shown in Ref. |30j (see also the Appendix of Ref. instead of finding A (r), T A (r), 
one can use the redundant variables b^(r), S^ u (t), with the following Dirac brackets assuring 
the validity of the orthonormality condition r] fJi1 ' — b^r] Ab b u ^ = [C^ al3 = TLv^rf"^ +r]^r]gr] L ' a — 
V"VsV fia ~ V^VsV" 13 are the structure constants of the Lorentz group] 

{sr,b p A y = v p %-v p < 1 b u A 

{ST, Sf}* = C^ a/3 S] 5 , (57) 
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so that while 7i^{r) pa has zero Dirac bracket with itself and with 7i^ u {r) pa these last 
six constraints have the Dirac brackets 

{7iT (t), n a ^r)}* = G^W\t) pa 0. (58) 

Let us now restrict ourselves to configurations with p 2 s > and let us boost at rest with 
the Wigner boost L fM u (p° s ,p s ) the variables 6^, Sf of the non-Darboux basis 

T M „M QV» n? K f 

•^si fsi u Ai u s i 'li i % 

of the Dirac brackets {.,.}* [in Refs. this step is missing; the final results in those papers 
are not changed, but the reduction associated with the following Eqs. (|62|) is not sufficient 
to get them as erroneously stated in Refs. ||]. The following new non-Darboux basis is 
obtained {x^ is no more a fourvector) 

~M M _i_ A/ ( \\ 9e p( U (Ps)) QUP 



^ r= ^\ps»S7 + v^p 2 s (s°/-s s , n 

V\JP 2 s{P°s+ r l\Jpf) 

1 - a S Ar v r P'i + 2?7 s a/p?77 mo - 

< -j=[Va{S° s s P ° ) + ' V ' S° s r pl 

VsJp 2 s P° + VsJp 2 s VsJp 2 (P°s+VsJp 2 s) 



p^pt 



V f 

r f 

K i — K i 



= sir - r>(p^ AB ( ^^^ p: - de *^ s)) P ^ = 

= sir + —j==7^ — rr>^Ps - sftf) + v^stp: - sfpff)] 
vyptKPs + v^p;) 
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Jr = - 5X + , (59) 
We have [cf. Eq.Q] 

{§l j , b T A }* = (6 ir 6 js - 5 ls 5 jr )b s A 

{sf , sfy = c^s] s , (60) 

and we can define 

N 



S? B = e A Mp s ))e B v H Ps ))S^ « [6?(r)6f - 6? (r)6?] E^iWW'"? + ^( r ) + 

i=l 

A 7 " 

+ [tf (r)&f (r) - bf(r)bf(r)} £^(t)/<(t). (61) 



i=i 



Let us now add six more gauge-fixings by selecting the special family of spacelike hy- 
perplanes orthogonal to p^ (this is possible for p\ > 0), which can be called the 'Wigner 
foliation' of Minkowski spacetime. This can be done by requiring (only six conditions are 
independent) 

t a(j) = b W - ^U(«(p.)) * o 

=► 6j(r) = 4«PsW a {t) « ^. (62) 

Now the tetrad 6^ has become e^(w(p a )) and the indices 'f' are forced to coincide with 
the Wigner spin-1 indices 'r', while o = r is a Lorentz-scalar index. One has 

N 



8=1 

JV 

+ (^-^V)E^)«?W 



i=i 



jV 

i=i 
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N I 

i=i 



rir crs s 

J? » xV s - £X ^h- (63) 

P° s +Vs\/Pt 

The time constancy of T| ~ with respect to the Dirac Hamiltonian of Eq. fl54|) gives 

^(r) - <(«(p.))] = mr) - e?( W (p s )),F D }* = 

= l~X aP {r){b^r),S sal3 (T)Y = A^(r)6, a (r) « 

A^(r) « 0, (64) 

so that the independent gauge- fixings contained in Eqs.(|62D and the constraints , W w {f) ps 
form six pairs of second class constraints. 

Besides Eqs. flHD , now we have [remember that x%(t) = — A M (r)] 

p = K = ^(Ps) 



<(r) = x?(t) = vG0V*(P«) -i w (T)eJf(tt(p.))<(«(p.)) 
</(t) = [s^KCp.))] 2 

SW = is M ^ s e^(ti(p s )), # rs = -5 rs H 7-r-i aM e(J(u(p a ))ar ai/ e^(u(p 8 )). (65) 

On the hyperplane S^(t) all the degrees of freedom z^(t,<t) are reduced to the four 
degrees of freedom x%(t), which replace x^. The Dirac Hamiltonian is now He = A M (r)7Y M (r) 
with 

N N 



H»{t) = ft - ti"(p.) £ r^m? + W - e?(u(p.)) E « 0. (66) 
t=i i=i 

To find the new Dirac brackets, one needs to evaluate the matrix of the old Dirac brackets 

of the second class constraints (without extracting the independent ones) 
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c 



{H a P,W s Y « 



{ft Q/3 ,T|}* = 

= SB B b aP e%HVs)) - V aa 4(<Ps))} 
{T P A ,W S Y = {Tj,T|}*=0 

{ = Sju[rr^iP»))-^l{u(p a ))] 



\ 



(67) 



Since the constraints are redundant, this matrix has the following left and right null 

f- ■ \ ( o ^ 



eigenvectors: 



0>af3 — apa 




[a a p arbitrary], 



c *Hps)) 



. Therefore, according to Ref. 



33| , one has to find a left and right quasi-inverse C, CC = CC = D, such that C and D 

\ 
J 



have the same left and right null eigenvectors. One finds 

/ " liVyr^f (u(p s )) -V6rtf(u(p s ))] 



C 



CC = CC = D 



Qaf3D 



Q BD 



\ 



i(^-e D "(n(p s ))e^ T ( M (p s )) 



(68) 



and the new Dirac brackets are 



1, 



{A,B}** = {A, B}* - -[{A,W s }*[r] yT ef(u(p s ))-T] ST e^(u(j) s ))}{T^,B}* + 

+ {A,TZ}*[r] au e*(u(p s )) - ^e* (u(p s ))]{H^ , B}*). (69) 

While the check of {7i al3 , B}** = is immediate, we must use the relation b^T^e Dp = —T^ 
[at this level we have T\ = T%\ to check {TjJ, B}** = 0. 

Then, we find the following brackets for the remaining variables x^,p^,7]\, k\ 

and the following form of the Poincare generators 

P^s 

Jg = XgP s — x s Pg + Sg 



(70) 



S° 



$13 = prp'S?. 



(71) 
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Therefore, x^ is not a fourvector and ffi, Hi transform as Wigner spin-1 3-vectors. 

Let us come back to the four first class constraints W l (t) « 0, {7tI m , 7i u }** = 0, of 
Eq. (|66D . They can be rewritten in the following form 



N 



N 



8=1 



H{t) = u^(p s )Hh(t) = rjs\Jp 2 s - ^rax/in: 

N 







(72) 



i=l 



The first one gives the mass spectrum of the isolated system, while the other three say 
that the total 3- momentum of the N particles on the hyperplane S^(r) vanishes. The 



Dirac Hamiltonian is now H D = X(t)H(t) — A(r) • 7Y p (r) and we have x s = {x^,Ho}** = 
— A(r)w M (p s ). Therefore, while the old x^ had a velocity x^ not parallel to the normal 
P = u^(p s ) to the hyperplane as shown by Eqs. (|65|) , the new x^ has x^ and no classical 
zitterbewegung. Moreover, we have that T s = I ■ x s = I ■ x s is the Lorentz-invariant rest 
frame time. 

Let us do the following canonical transformation [with the j a i of Eqs.(|6]), (|7|)] 





Ps" 


Vi 









T 

- 1 s 


z s 












-4 






Pa 


K 




V+ 







(73) 



rp Ps ' Ps ' 3^8 

Vs\Jp1 Vs\JpI 
e s = Vs\ff s 

-> r~oi ~ P s ~o\ 

Zs = Vs\/Pt{X s ~ —X s ) 
Vs\P 2 s 



N 



v+ 



N 



8=1 
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N 
i=l 



N 



-I 



Pa = VNj2^iVi 
i=l 

1 N 

7T„ = 



whose inverse is 



x s 

P°s 
Ps 



Nfi 



e., e s 



iV-1 



Vi = V+ + 
1 



0=1 

AT-1 



By using the constraints k 

N-l 2 
o=l 

AT 



N 

0, we get 



0=1 



o=l 



s: 



i=i 



\ 



N-l 



N-l 



N 



Tin 



a=l 



a=l i=l 



(74) 



(75) 



iV-l 



irr, 



+ N(J2%^'a) ■ (76) 



a=l 



Note that we have N+l 3-vectors z s ,ff + ,p a for the N particles, but that the constraints 
k + « imply that the ff+'s are gauge variables. We can decouple these 3 canonical pairs 
from the others: the simplest way is to add the gauge-fixings ff + « [let us remark that 
z M (r, 0) = x^(r), i.e. this is the origin of the coordinate system on the hyperplane from 
the point of view of Minkowski spacetime]; so we remain with the correct number of N-l 
relative variables p a ,Tr' a plus the center-of-mass ones x^,p^ and with [compare with S or = 
- Ell 1 p r a e Ra at T Ra = 0, eq.fllj)] 
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1 N-l N 



JV-1 ' AT-1 

2 



+^(E7«<) =-Epl r ^- ( 77 ) 



a=l a=l 



where we anticipated a result of next Section. 

In this way we obtain the rest-frame instant form T s = const, of the reduced problem 
for the N-body system with the center-of-mass motion separated out. This constant motion 
is taken into account by the observer which looks at the hyperplane Y,w(t). The evolution 
in T s of the reduced system is governed by the remaining constraint 



N N-l z 

W(r)=e.-5> rn? + JV(£7-i<) = C - H™ « 0. (78) 

i=l \ a=l 

If we would add the gauge-fixing T s — r m and would go to Dirac brackets, we would 
obtain a Hamilton- Jacobi description in terms of Jacobi data. If we want to reintroduce 
an evolution for the relative variables, the natural time parameter is T s and the associated 

(T ) 

Hamiltonian is e s = H R ; the Jacobi data would then become the Cauchy data for this evo- 
lution. In the nonrelativistic limit Hj[ would go (modulo m^c 2 terms) in the Hamiltonian 

Hr of the reduced N-l relative variables problem: H = =%j h Hr would be the total 

Newtonian Hamiltonian. 

The vectors z s , k s describe the frame-dependent position of the canonical center of mass; 
in the rest frame, P % = e s (l;Q),z s defines a point in the hyperplane Y> w (r) orthogonal 

o 

to P which can be used to build the covariant noncanonical Fokker center of inertia. 
Therefore, the worldtube defined by the frame-dependent canonical center-of-mass positions 
( s = Zs/rjsJp 1 * arises naturally also in this 1-time descritpion (the rest-frame instant form 
T s = const.). Let us note that, while e s = n s ^Jp^ is determined by the invariant mass of 
the physical system on the hypersurface [see Eq.([7|)], k s describes the orientation of the 
hyperplane with normal f = u^(p s ) embedded in Minkowski spacetime with respect to an 
arbitrary observer: k s is the collective velocity of the physical system as seen from outside 
the hyperplane and not the three-momentum, determined by the energy-momentum tensor, 
on the hypersurface, which is the vanishing k + = 0, since the hyperplane corresponds to the 
rest frame. 
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For generic masses rrii, Eq.flTBD describes all the 2 N branches of the mass spectrum by 
considering all the possible combinations of rji = signp° = ±, i.e. the solutions in e of \ ~ 
of Eq.fl26|) when 6R a w 0. In the 1-time theory there is a different action ( |38| ) associated 
with each branch of the mass spectrum of the N-time theory. 
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IV. THE CONNECTION BETWEEN THE 1- AND N-TIME DESCRIPTIONS 



We have arrived at two descriptions of N scalar free particles: i) the 1-time theory 
with variables x^, p a , 7r , ff + , k + and with the first class constraints (|72|), i.e. k + ~ 



and H = e s — J2f=i Vi\f m i + N(Y<a=i lai^' a ) 2 ~ 0; ii) the N-time theory with variables 
x 1 , p a , 7r a , Tna and with the first class constraints eR a ~ plus a mass-spectrum constraint 
X ~ 0, Eg. . The last line of Eqs.(PT|), giving the original constraints fa ~ 0, shows that 
it is the mass-spectrum constraint x ~ that governs the dynamics in both descriptions, 
since 7i ~ identifies one of the 2 N branches of x ~ 0- 

To compare the two descriptions, we have to identify two isomorphic phase spaces: i) in 
the 1-time theory, we add the gauge-fixings fj + « [natural from the comparison of Eqs.(|l5|) 
and (|77|) 1 , so that the first phase space is spanned by x 1 ^ , p^ , p a , n a , with the constraint 7i w 0; 



ii) in the N-time theory we add the gauge- fixings Tfj a « implying Eqs. (|33|) [so that the 
final variables of the second phase space are x^,p^, p a , n a ], and we choose among the 2 N 
branches of the mass spectrum the one corresponding to the given choice of the signs rji of 
the 1-time theory. 

Since we do not know the expression of the original variables xf,pf in terms of the 
variables of Eqs.([H]), to find the relation among the variables of the two phase spaces we 
have initially to go back to Eqs.(|9]) 

1 N-l -I N-l 

xt = x» + -= £ laiK =^+-7hE %A(<p))e*(u(p))R: = 
VN a=1 VN a=1 

1 N ^ v» 
= x» + -= Y: l^{-T Ra + e?(u(p))p r a ) -> 

j N-l 

^7^=^=021^=^=0 + -7=e?(u(p)) TaiPa 

V^V a =l 

i N-l 

= -t-jp" + vn £ %m = 

ly a=l 

AT-1 AT-1 

= (^7 + v / iVE 7«iCitoK(p) + Vfa£(«(p)) £ 7ai< - 

iV a=l a=l 

. . ? 



(79) 
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On the other hand, Eqs.(|39D, d62|) and ( |66|) imply the following identification with the 
variables on the hyperplane orthogonal to p^ 



1 de B (u(v )) 1 Ar_1 

* °Psfi VJV a =l 

= K + — f^T~ r=T IP-^ + - S °s U ^)] + 

V\/Ps(P°s+V^P 2 s) P s 
i JV-1 
+ -7=e^(w(p s )) £ 7a,Pa = 

viv a=1 

Vs y/Pt Ps + Vs y/pi Vs ypj {P°s + Vs\/Pi) 
1 JV ~ 1 



Pi\T Ra =e Ra =o = Vi\/mf + k^(t)u^(p s ) + e£(u(p s ))<(r) 



JV-1 



^mj + ^(r)u^(p s ) + e%(u{p s ))[—K f + + VNY, 

iV a=l 



|k+=0 



JV-1 



N-l 



in: 



+ 7-<W) u"(p.) + Vfo(f(u(p.)) £ 7^. 



a=l 



a=l 



We see that |T fla =e fla =o satisfies automatically pf = mf and that 



N 



N 



1=1 



8=1 



1 



JV-1 



a=l 



if we have 



(80) 



i + N(Y: 7«<(r)) ^ + <£(«(?.))«; « (81) 



A? 



i=l 



\ 



JV-1 



m: 



+ JV(E»)) 



a=l 



by using Eqs.(|59|). Moreover, we have 



(82) 



JV 



Qa\T Ra =i Ra =0 — —7ff^2laiPi\T Ra =£ Ra =0 



N 



1=1 

N 



N-l 



nr, 



+ N(Y: 7«i<(r)) «"(p), 



0=1 



(83) 
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so that for T Ra = e Ra = k + = we find 



1 N A 

e Ra\T Ra =e Ra =0 ~ ~7= X! 7aiVi 



\ 



iV-1 



m : 



(84) 



a=l 



Equation fl34|), with Eq.(|8|), implies that the second one of Eqs. (|79| ) becomes the second 
one of Eq.(|^) with 7?^ = n a . Equation fl34|) gives all the looked for solutions of Eq.(p6[) 
[or of Eqs. (37)] for T Ra = e^ a = according to which branch of e ~ e s has been chosen in 
Eq. fl32"|) , as can be explicitly checked. 

Moreover, for each branch the mass spectrum, given by Eq.(^) for a given choice of the 
signs rji, we have the following solution of Eqs. (|3"TD without the restriction T Ra = e Ra = 
[remember that 7? a = 7? a = n' a ] 



N 



J2 IfaiVi 



i=l 



\ 



N-l z 
a=l 



2 1 ' 



(85) 



if Eq. (^) is rewritten in the weakly equivalent form 



v 

e ~ e s ~ XI ^ 



\ 



N-l 



2e 



N-l 



m l + N ( 7<m<(t)) + -7= Yl IK*™- 



a=l 



(86) 



b=l 



On the branch defined by Eq. (|86|) , Eqs. (|85|) imply the following form for pf of Eqs. (ffiH) 



Pi « 



\ 



iV-l 



iV-1 



2e JV '~ 1 — 
2 + # ( E 7«i<(*r)) +^=E 7^i» + v 7 ^ X 7a*P • vr a ] 



AT-l 



a=l V A' 

AT-l 



6=1 



a=l 



ft 



e p° + e 



a=l 



N-l 



2e 



7V-1 



and Eqs. (|26|) give 



ml + N( E 7ai<(r)) + -= E 7wej»] 



a=l 



(87) 



b=l 



2e iV " 1 

E 7ai6i?a ~ 0, 



(88) 



a=l 



consistently with \ ~ due to Eq. (|8~6"D . Therefore, one has [n^ means a product over all 
the 2^ choices of the signs rji] 
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N N-l 

2 



Vi 1=1 



i=l a=l 
N 



„ + iV(E 7«<(r)) +^=E fa*J»] « 0. (89) 

\ a=l V 6=1 



Even with Equations (p5|), (pq) , we cannot get the explicit inverse canonical transfor- 
mation giving of Eqs.fl20D in terms of ,p^,T Ra , e Ra , p a ,7t a = Tx a of Eqs. (j3ll) inde- 
pendently from the choice of the branch of the mass spectrum, since Eqs.(|85|) are a branch- 
dependent solution of Eqs.(j37j). 

While Equations (|87D give the answer for the second half of Eqs. (|79|) , the discussion of 
the first half is more involved. 

By using Eqs. (pi]), fl33|), (p8|), ( |i5|) evaluated at T Ra = e Ra = 0, we get 



e(p° + e) e(p° + e) 

1 N-l 1 



a=l 6 



1 ~ n^S sr v r 



so that Eqs.(^) and (|80|) imply at T^a = e/?a = 



iV-l 



< « a" + -7=e?(«(p)) E 7a*P a 

VJV a =l 



+ -<(«(p))5* - + ^(t*(p)) E %iP r a 



By comparing the two expressions for /x = and for /j = k, we find 



(91) 



Pa = Pa (92) 
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due to Eqs.(|76f) and ( fT7|) [for p a = p a we have S° r = S or and Sl s = S rs ], consistently with 
the identification of the two reduced phase spaces, whose final variables can be chosen as 
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V. THE MASS SPECTRUM AND THE INTRODUCTION OF INTERACTIONS 



In the 1-time theory, by varying r\i = signp° in Eq.([78|) we obtain the 2 N branches of the 
mass spectrum of N free scalar particles in the rest-frame instant form, when p 2 > and the 
masses m^'s are generic, so that degeneracies are avoided. For each of these branches there 
is a well defined Lagrangian, which gives rise to the constraints of Eqs.(^) on arbitrary 
spacelike hyperplanes: only the configurations with p 2 > can be reduced to the rest- 
frame instant form. Being in classical relativistic mechanics, we must take into account also 
negative masses and energies. 

To see which other configurations may be allowed, we have to revert to a study of the 
N-time theory with its first class constraints 4>i = p 2 — m 2 « 0. 

Let us consider in detail the case N=2; when p 2 > and rji = +, i=l,2, we are considering, 
in absence of interactions, the special configuration of the kinematics of forward s-channel 
eleastic scattering, in which we have s = (pi + p 2 ) 2 = p 2 = e 2 , t=0, u = (pi — P2) 2 = ^Q 2 as 



special values of the Mandelstam variables [see Ref. ||32|| ], so that, given the following two 
combinations of the fa's [see Appendix B] 

X+ = 2(0i + <j> 2 ) =p 2 + AQ 2 - 2{m\ + ml) « 

X- = \{<f>x-<h)=V'Q-\{m\- m 2 2 ) w 0, (93) 

we find that x+ ~ is equivalent to the fundamental relation s + t + u = 2{m\ + m|). 
Let us discuss the various possibilities of choice of the value of the masses. 



1) m\ 7^ m 2 (mi > m 2 > 0). In this case, the first class constraints (|93| ) define a 
constraint submanifold 7 (coisotropically embedded in phase space [|]]), which is the disjoint 
union of three strata with p 2 > (the main stratum dense in 7), p 2 = and p 2 < 
respectively [the stratum with p M = is excluded by X- ~ 0] . 

la) Stratum with p 2 > 0. According to Appendix B, one has [Q 2 ~ Q\ = —tt 2 < 0; 



2 2 
l- m 2 I 



Q 2 - Q 2 , 



X- = P ■ Q w 
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x+ = + 4Q 2 ± + 4^j^ - 2(m 2 + m 2 2 ) ^ p 2 + AQ 2 + ^ - m ' - 2(m 2 + mj{) = 
= ^(p 2 -M 2 )(p 2 -M!)^0. (94) 



with M-t given in Eqs.([I 

There are four branches for e = rjy/p 1 in terms of 7? 2 = — Q\ pa —Q 2 



e pa ±{\fm{ + 7T 2 + ymi + vr 2 ), ^1=^2 



e « ±| (ymf + 7f 2 - \jml + 7F 2 , | r/i = -rj 2 (95) 

For 7T 2 — > [corresponding to <5 M = (pj* / / P2); Q 2 = is excluded by X- ~ 0] we have the 
(s-channel) thresholds e pa ±(mi + ^2) and (u-channel) pseudo-thresholds e ~ ±| mi — m2 |. 

lb) Stratum with p 2 = 0. It exists for r\\ = —i]2: for pf timelike, i=l,2, we have 
= ( Pl + P2 ) 2 = (771I +772I P2I) 2 — (Pi +P2) 2 ~ {viy/™%. +Pi + ^2\/ m 2 -(pi + P2) 2 = 



m 2 + m 2 , + 2rjir]2\fm 2 + p!ym 2 , +j% - 2pi ■ p 2 p> 2=0 - 277 1 77 2 7Ti 2 y + pf = m 2 + m 2 , 
f?if?2 = -• 

Now the constraint % + « implies Q 2 ~ |(m 2 + m|) > and two branches 



Q° ^±Ji(m 2 + m 2 ) + g 2 (96) 



of excitations of relative energy with total energy given by 

1 r l 



P^Q-J^rnl-mD+p-Q]. (97) 



lc) Stratum with p 2 < 0. It exists for r)\ = —1]2', the constraint x+ ~ implies Q 2 ~ 
^[2(m 2 + 7T1 2 .) — p 2 ] > 0. In terms of = Q 2 — ^2 ~ Q 2 — ( ' r " 1 4p T' 2 ' > , we have two tachionic 
branches 



P 2 «-(V^i-^±V^i-^) ( 98 ) 

2) m\ = 1712 = m > 0. We have 

X+ = p 2 + 4Q 2 - 4m 2 pa 

x _ = p . Q « 0, (99) 
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and there are three strata with p 2 > 0, p M = 0, p 2 < respectively [p 2 = is excluded by 
Q 2 w m 2 > and p ■ Q ~ 0]. 

2a) Stratum with p 2 > 0. We have Q 2 ~ Qi = — v? 2 < and there are two branches 

e « ±2Vm 2 + 7T 2 , [e « ±2mforpf = p%, i.e.Q* = 0], (100) 

which are the limit for mi — m 2 — > of the two branches with rji = r] 2 of the p 2 > stratum 
with mi 7^ m 2 . 

2b) Stratum with p^ = [pi = — pg =>- r\\ — —r] 2 ]. Since we have the constraint 
X+ = 4(Q 2 — m 2 ) « 0, there are two branches 

g° « i^ + g 2 (ioi) 

of excitations of relative energy: they correspond to two degenerate e = levels obtained 
from the limit mi — m 2 — * with r\\ = —772 of the p 2 > stratum with mi 7^ m 2 [see Ref. 
p4| for the spin 1/2 case]. 

2c) Stratum with p 2 < 0. It exists for rj\ = — r/ 2 [p 2 < is equivalent for p 2 = to 
2m(m + r]ir]2\Jm 2 + pf) < 0], and we have Q 2 — |p 2 > 0, so that there is a tachionic branch 

p 2 « -4(Q 2 -m 2 ) « -4(Q 2 L -m 2 ). (102) 

3) mi = m 2 = 0. The constraints p 2 ~ 0,p 2 ~ give 

X+=p 2 + 4Q 2 ^0 

X-=p-Q~0, (103) 

so that there are three strata with p 2 > 0, p 2 = 0, p M = 0, respectively [p 2 < is excluded, 
because it would imply p 2 = (?7i|pi| + 772IP2I) 2 - (pi +P2) 2 = 2\ Pi\\ P2KV1V2 - cos6 12 ) < 0]. 
3a) Stratum with p 2 > 0. It has Q 2 ~ = — v? 2 < and there are two branches (with 

m = m) 

e « ±2^ (104) 
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tangent at vf 2 = (Q» = 0). 

3b) Stratum with p 2 = 0. It implies 2| pi\ | 1 ( 1 7i ? 72 ~ cos 6i 2 ) = 0, i.e. p\j jp 2 with r/i = r\ 2 
for #i2 = and with r/i = — r\ 2 for Q\ 2 = it. We have Q 2 ~ 0, namely two branches 

Q° « iv^ 2 " (105) 

of excitations of relative energy with total energy 

f » ^ „ ± ^|. (106) 

3c) Stratum with p M = [pi = —p 2 , i]i = —i] 2 ]. It has Q 2 ~ 0: there are two branches 
of excitations of relative energy 

Q° « iV^ 2 . (107) 

For N > 2 the discussion is more involved, since there are many degenerate cases with 
complicated discussions of the allowed Poincare orbits. 

For N=3, p 2 > and mi,m 2 ,m 3 generic, there are 8 branches 

e ~ ±[+\fm\ + 3(7n7?i + 7 2 ivf 2 ) 2 + \/mi + 3(712^1 + 722 v?2) 2 + 

+ \Jml + 3(7i 3 7Ti + 723V?2) 2 ] 

e ~ ±| + \/mf + 3(711??! + 721 vf 2 ) 2 + \jm 2 2 + 3(7i 2 7?i + ^ 22 n 2 ) 2 - 

- \/ m 3 + 3(7i 3 7fi + 723^2) 2 I 
e ~ ±| - \Jm\ + 3(7n7?i + 721 vr 2 ) 2 + \fm% + 3(7i 2 vfi + % 2 n 2 ) 2 + 

+ \/ m 3 + 3(7l37Tl + 7237?2) 2 | 



e w ±| + \fm 2 + 3(7n7?i + 7 2 i7f 2 ) 2 - \Jm 2 2 + 3(712^1 + 7227T2) 2 + 
+ + 3(7i 3 7?i + 72 3 vr 2 ) 2 | . (108) 

In the limit mi = m 2 = = m, the first two branches start from e ~ ±3m, while the other 
three positive (negative) branches are tangent at e « m (e w — m). 
This pattern is valid for all N=2k+1. 
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For N=4, p 2 > and mi,m 2 , m^, ""U generic, there are 16 branches 
e ~ ±[+y/m% + 4(7n7?i + 7 2 i7? 2 + 7 3 i7?3) 2 + + 4 (7i27fi + 722*2 + 732^3) 



2 + 



+ v/m| + 4(7x37?! + 7 23 7r 2 + 73 3 7r 3 ) + J ml + 4(7 14 7Ti + 7 24 7f 2 + 734^3) 



e ~ ±| + \Jm\ + 4(7n7fi + 7 2 i7r 2 + 7 3 i7r 3 ) 2 + \[m\+ 4(7i 2 7Ti + 7 22 7? 2 + 7 32 7f 3 ) 2 + 



+ V m i + 4(7i3vri + 723^2 + 733^3) - V ml + 4(7i 4 7Ti + 7 24 7r 2 + 734^3) 



e ~ ±| - \Jm\ + 4(7n7Ti + 7 2 i7r 2 + 73i7r 3 ) 2 + \fm% + 4(7i 2 7Ti + 7 22 7r 2 + 7 32 7f 3 ) 2 + 



+ V m 3 + 4(713^1 + 723V?2 + 7337?3) + V m 4 + 4(714^1 + 724^2 + 7347?3) 



±| + \Jm\ + 4(7n7fi + 72i7f 2 + 7 3 i7f 3 ) 2 - \Jml + 4(7i 2 7Ti + 7 22 7r 2 + 7 32 7f 3 ) 2 + 



+ V m\ + 4(7i 3 tTi + 7 23 7r 2 + 733 ^3) 2 + \j ml + 4(7i 4 7Ti + 7 24 7r 2 + T^^s) 2 I 
~ ±| + \jm\ + 4(7n7fi + 7 2 i7f 2 + 7317T3) 2 + ^/m| + 4(7i 2 7fi + 722*2 + 732^3) 
- \jm\ + 4(7137?! + 7 23 7r 2 + 7337T3) 2 + \/ml + 4(7i4*\ + 724*2 + 734 7?3) 2 I 
~ ±| + \fm 2 + 4(7n7?i + 7 2 i7? 2 + 7 3 i7f 3 ) 2 + \Jml + 4(7 12 7Ti + 7 22 7r 2 + 7 32 7r 3 ) 



'ml + 4(7i 3 7fi + 7 23 7f 2 + 733VT3) 2 - V m 4 + 4(7l47Ti + 7 24 7T 2 + 734??3) 2 | 

e ~ ±1 - + 4 (7n^i + 72i*2 + 73i7?3) 2 + \jm\ + 4(7x27?! + 7 22 vf 2 + 7 32 7r 3 ) 2 + 



+ \jm\ + 4(7i 3 7Ti + 7 23 7f 2 + 733^3) - \]ml + 4(7i 4 7Ti + 7 24 7r 2 + 734^3) | 
e ~ ±| + \Jm\ + 4(7n7Ti + 7 2 i7? 2 + 7 3 i7f 3 ) 2 - \Jml + 4(7i 2 7Ti + 7 22 7? 2 + 7 32 7f 3 ) 2 + 



+ V m 3 + 4(713^1 + 7237T2 + 733^3) ~ V m 4 + 4(714^1 + 7247?2 + 734^3 ) | . (109) 

In the equal mass limit, there are two branches starting at e ps ±4m, four positive (negative) 
branches starting at e ps 2m (e ~ —2m) and six (three positive and three negative) branches 
starting from e ~ 0. 

This is the pattern for N=2k. 

As we have seen, already in the free case there are generically tachionic strata corre- 
sponding to spacelike Poincare orbits. At least at the classical level, these strata have to be 
excluded to avoid problems with Einstein causality. Instead, the rest-frame 1-time instant 
form is, by construction, free of these tachionic strata. 

45 



For N=2, strata with p 2 < appear always when 771 = —772; this fact will extend to 
N > 2. Now there is the problem of the interpretation of particles of negative mass and 
energy (rji = — ); they are present in classical relativistic mechanics due to the two branches 
of timelike Poincare orbits and have a well defined nonrelativistic limit, connected with the 
negative mass representations of the extended Galileo group 



If we make the first quantization of a scalar particle, the first class constraint p 2 —m 2 m 
goes into the one-particle Klein-Gordon equation (□ + m 2 )(j)(x) = with <fi(x) a complex 
wave function. It is known |36| (see also Ref. [pCfl ) that this equation admits two kinds of 
scalar products: i) the standard nondefinite positive one (4>a,4>b)i = / d^a^ 4>\(x)^ 
[d^ =d Ai — d M ] associated with the conserved current ji(x) = |0*(x)<9 M 0(x), which is inter- 
preted as a charge (electric charge, strangeness,..) current with the sign of the charge corre- 
sponding to positive- and negative- norm (i.e. energy) states; ii) a positive definite nonlocal 



one [<pA, 4>b)2 = I d 3 cTfj, j<f>A( x ){l ~ [77 = id° / ym 2 — d 2 is the nonlocal Lorentz- 

scalar (under proper Lorentz transformations) "sign of the energy" operator satisfying 
[rj,x ,M ] = [fj, id^ 1 } = 0], associated with the conserved current j% (x) = z<f>*(x)(f) — fj^d^^x). 
The negative mass (and energy) states are interpreted as describing antiparticles: for charged 
scalar particles (7r + 7r~, K + K~, K°K°) positive- (negative-) energy states describe particles 
with charge +Q (antiparticles with charge —Q), which can be connected by a "charge con- 
jugation" operation <p{x) 1— > i] c <f)*(x), \r] c \ = 1, under which the action is invariant [for the 
electromagnetic potential one has A^x) 1— > —A^x)]; for neutral particles like n°, coinciding 
with the antiparticle, see Ref. |37 . 



In Ref. [Q, it was shown at the classical level that, since the proper time of a particle 



is defined by ds = ±Jdx fJ "q^ u dx u if dx° < 0, we have ds = ±mdT J x^r) x u , so that dr > 



corresponds to dx° < if x° = ^ < 0; therefore, an evolution towards the future in dr for 
a particle with negative mass -m corresponds to an evolution towards the past in dx° and 
this is in accord with the use of the complex Stueckelberg- Feynman Green function Gf(x) 
[see for instance Ref. P9| ], which propagates the positive- (negative-) frequencies forward 
(backward) in x°, in association with the first-quantized Klein-Gordon equation. Finally, the 
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equations of motion of a particle of mass m and electric charge e in presence of gravitational 
and electromagnetic fields depend only on the ratio e/m, so that an antiparticle of mass -m 
(propagating forward in r but backwards in x°) and charge -e satisfies the same equations 
as the particle (m,e) if parametrized in r, in accord with the charge conjugation invariance. 

This is what happens with the parametrization in r of the 1-time theory. Let us remark 
that in the rest-frame instant form one has {e = rjy/p 1 ^ = p '^— } = 1, so that rjm is 

r\\l p 2 



associated with 77 1 T\. 

In contrast, in the Klein-Gordon quantum field theory, a hermitean quantum field 4>(x) 
[quantization of a classical real Klein-Gordon field 4>(x) = <fi*(x) for which the complex 
Green function Gf{x) cannot be used, but only the real retarded and advanced ones] has the 
positive- (negative-) frequency part associated with the creation (annihilation) operators of 
a scalar particle with positive energy For a complex quantum field 4>(x), ^(x) [quantization 
of two real Klein-Gordon fields, 4>{x) = -^(<pi(x) + i(p2 (%))} the two kinds of creation and 
annihilation operators are associated with a particle of mass m and charge +1 and with an 
antiparticle of mass m and charge -1 respectively. 

The conclusion of this discussion is that, in the classical background of relativistic particle 
physics, we have to consider only the branch of the mass spectrum of N scalar particles 
with all the masses positive and with particle and antiparticle distinguished by opposite 
charges [only neutral scalar particles cannot be described in this way, but they (by the 
way also the charged scalar ones) are supposed to be bound states of spin 1/2 quarks]; the 
pseudothresholds of the lower positive branches are connected to the thresholds of the other 
existing kinematical invariants (relevant for scattering theory due to the crossing property). 

Let us now consider the introduction of action-at-a-distance interactions. This problem 
has been studied in the N-time theory (see the bibliography of Refs. [[40|--f42"|| ) . In the frame- 



work of models with N first class constraints, a closed form of the constraints is known for 
N=2: this is the DrozVincent-Todorov-Komar model [pi|-f45|| , on which Ref. |2(J is based. 



The constraints take the general form (only the case V(R\) has been studied in detail; see 
Refs. |2(|,[l6j for the nonrelativistic limit) 
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p *-m? i+ V(Rl,p 2 ,Q 2 ,Q-R x )nO, 

{01,0 2 } = O. 



1,2, 



110) 



By construction we have p 2 > 0. The potential V may be either confining or separable 
(see Ref. [|(J for a finite-range potential, in which the particles see each other only if their 
Cauchy data are restricted to be compatible with the first class constraints in the interacting 
region) . 

For N > 2 there are models with a closed form of the constraints (pi ~ 0, {(pi, (pj} = 0, 
only for confining potentials (see for instance Ref. [19|]). To find separable potentials such 
that {(pi, (pj} = requires the solution of complicated nonlinear partial differential equations 
connected with the necessity of introducing 3-, 4-,.. N-body forces pTFfl , and no one succeeded 
in solving them; it was only shown that solutions exist, in which the first class constraints 



are expressed as series in the coupling constants []48f and that the solution is unique [49 



See also Refs. [p^Jlql - 

In contrast, in the 1-time theory (the rest-frame instant form for p 2 > 0), we can intro- 
duce the interactions like in Newtonian mechanics and use the nonrelativistic definition of 
separability of the interactions. It is the transition from the 1-time to the N-time theory 
that contains all the previous difficulties. 

The most general form of the final constraint (|78|) of the 1-time theory is 



H 



N 

52 m 

i=l 



\ 



N-l 



1..N 



m\ + Vi + N(Y, larfa) ~ 52 U H = C * 



H 



(Ts) 
R 







a=l 



N-l 

V i = V A52 (^ah - lak)pa] 
a=l 
N-l 

a=l 



111) 



With Vi = and N=2, it is the form prescribed in Ref. [|50| [see also the bibliography of 
Ref. pl |, but, as we shall see in the next Section, the 2-body Coulomb potential coming from 
the longitudinal modes of the electromagnetic field is of the type U^; therefore the models 
with U-type potentials can be thought as deriving from couplings to gauge field theories. 
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The action-at-a-distance interactions (instantaneous in the rest frame) of models like the 
one of Refs. [^3|-p5| (for instance the relativistic harmonic oscillator) are of the type (addi- 
tive to mj). If Vt = Y^fi V ik [Y,^{%i-%k)M+Y} k ;^i V ikh [^=K%i-%k)Pa, E^ilai - 
lah)Pa\ + ..., we have 2-body potentials Vik, 3-body potentials Vn-h and so on. The potentials 
are separable if, whenever particle i decouples from the other N-l (the same must hold true 
for clusters of particles), we have V$ — > and C/y — > 0, so that no one of the surviving 
potentials depend on the index "i" . 

Let us remark that we can have two different 1-time theories with the same non- 
relativistic limit H R : i) = ££i + V t + N(Ea=x Vaf\ n) = 

Eili m VW + N(Ea=i %rfaf + U if U = E,=i instead, the N-time theory with a 
unique potential |1| \pj - mf + j^V « 0] would give N nonrelativistic constraints [see Ref. 
pj for the case N=2] ^ = E,-^--^ « with £^ = E-H = E- % ~ +H R « 
with the same -£/r of the 1-time theories if C/ = V\ 

Therefore, there are two 1-time models for the relativistic harmonic oscillator 
when N=2: i) the 1-time version of the DrozVincent-Todorov-Komar model H R = 
Viyjml + Vip*) + 2f 2 + r] 2 \jm\ + V^p 2 ) + 2vf 2 ; ii) a model with ^ Ts) = r/i^mf + 2t? 2 + 
r]2yrn% + 2tt 2 + U{p 2 ) with [/ = j\V [// = r ^ 1 _^ 2 is the reduced mass]. The second model 
should be interpreted as coming from a gauge field theory producing action-at-a-distance 
interparticle harmonic forces. 

Let us also note that the number of branches of the mass spectrum in the free case is 
a topological number: it is the dimension of the zeroth homotopy group of the constraint 
hypersurface, counting how many disjoint components are in it. While certain interactions 
preserve this number, generic interactions will change it; therefore, the interactions should 
be classified according to the (lacking) theory of intersections of noncompact hypersurfaces 
in phase space (when certain mass gaps disappear). 

In the approach of this paper, based on the theory of canonical realizations of the Poincare 
group in phase space, we cannot introduce external interactions, without destroying the 
technology we are using, which presupposes the existence of ten finite conserved Poincare 
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generators for every system under study. External interactions should be thought as limits 
from our isolated systems, when some its subsystem has the invariant mass tending to infin- 
ity; alternatively one can try to introduce them in the final reduced form of the constraints, 
like those in Eqs. fllll ). 
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VI. THE ELECTROMAGNETIC INTERACTION 



In Ref. ||21|| , the case of a charged scalar particle interacting with the electromagnetic 
field was considered. The 1-time theory was used because it is not known to us how to 
develop a covariant N-time description, in which the Poisson brackets of particle and field 
constraints can be evaluated in a covariant way. Here, the case of N charged scalar particles 
will be considered, with the electric charge of each particle described in a pseudo classical way 
by means of a pair of complex conjugate Grassmann variables [Q #i(r), 9*{t) satisfying 
[ij = I* = 9*9i is the generator of the U em {l) group of particle i] 

Q] = ef = 0, 9-fi* + 9*9i = 0, 

9 i 9 j = 9 j 9 i , 9 t 9*=9*9 u 9*9* = 9*9*, i^j. (112) 
This amounts to assume that the electric charges Qi = Ci9*9i are quatized with levels and 



ei [53|. 



On the hypersurface E(r), we describe the electromagnetic potential and field strength 
with Lorentz-scalar variables A a (t,<j) and F a ^(t,(t) respectively, defined by 

A A (r, a) = ^(t, a)An{z(T, a)) 

Fab^B) = d A A s (r,a)-d s A A (r,a) = zftr, <?)z%(t, a)F^(z(r, <?))■ (H3) 
The system is described by the action 



S = J dTd 3 aC{T,a) = J drL{r) 
L(t) = J d 3 aC{r,a) 



i N ■ 
i=i 

N 



^ 3 (^ - Vi( T ))h m i\/ 9tt(t, a) + 2g Tf (r, 0% f (j) + g ft (r, v)ri!( T )v!( r ) + 



+e l 9*{r)9 l {r){A T {r 1 3) + A f (r, 3)^(t))) - 
1 



[ v g(r,a)g^(T,a)g^(r,a)F AS (T,a)F CD (r,a), (114) 
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where the configuration variables are z^(r,a) A^(r,a), f)i(r), Oi(r) and 9*(r), i=l,..,N. We 
have 

-\V99 AC 9 6£) FabFcd = -sn\\^ F r?l f ~ SF T-s- T g 9Tvl if F rsl ™ 

The action is invariant under separate r- and cx-reparametrizations, since A t (t, a) trans- 
forms as a r-derivative; moreover, it is invariant under the odd phase transformations 
56, i i— > iaOi, generated by the Tj's. 

The canonical momenta are [E f = F fT and B f = e fS iF S t are the electric and magnetic 
fields respectively] 

oz T [t, a) i=1 

\J9tt(t, v) + 2g T f(r, B)t]1{t) + g f ~ s (r, a)f]t (r)r]f (r) 

+ ^^l(g TT Zr, + 9 Tf z f ;)(r, a)g A6 (r, a)g 6i) (r, a)F M (r, a)F CD (r, a) - 
- 2Mr, d){g' g'\r + g M 'g n g")i.. cf) + 
+ z f ,(r, a)(g Af g^ + g A ^g fd )(r, a)g^(r, a)]F M (r, a)F d£) (r, a)} = 



7r T (r, a) 
7r r (r, a) 



dL 

dd T A T (r, a) 
dL 



= 0, 



7(r, a) 



Y s (t, a)(E s (r, a) - g T ,(r, *) 7 w (r, 5)c«fS f (r, a)), 



K>ir(j) — 



dd T A f (r,a) ^(r, a) 
= i{t, dp 

dL(r) 



Y s (r,a)(F TS -g T ^F ilS )(T,a) = 



drfiij) 
= r/im 



g T r(r, ffi(T)) + g fS {r, ffi(T))rjf (r) 



V Wr(T,^(T)) + 2grf(T,Vi(T))flt(T) + Jf S (t, f]i (t) )?)[ (r)?)| (t 

+ e^*(r)^(r)A,(r,r/ i (r), 



+ 
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. , dL(r) i , . 
o9i{r) 2 

*»*W = 7^ = -ft( T )> (115) 

09* (r) 2 
and the following Poisson brackets are assumed 

{z»{T,o),p v {r,a'} = -^{a - a') 
{A A (T,a),n*(T,a')}=45 3 (a-a') 

{vi(T),K? j (T)} = 6 ij S n 

{9 i (T),n 0j (r)} = 

{e:(r),n e , J (r)} = -5 lJ . (116) 

The Grassmann momenta give origin to the second class constraints n ei + |0* w 0, 
77(9* i + PS [{7T6ij + |0*,7T0*j + = — 7^ and 7i>j are then eliminated with the 
Dirac brackets 

{A, B}* = {A, B} - i[{A, n e , t + ft}Ue* i + ft, B} + {A, 7r e * t + ft}{n ei + ft, B}] 

(117) 

so that the remaining Grassmann variables have the fundamental Dirac brackets [which we 
will still denote {., .} for the sake of simplicity] 

{9 l (r),9 J (r)} = {9*(r),9*(r)} = 

{9 l {T),9]{r)} = -i5 ir (118) 
Again, we obtain four primary constraints 
a) = p M (r, a) - Z M (r, a)[T TT (r, a) + 

Vi y/m1-Y B {r, SI(rj - e^*(r)^(r)A,(r, - e^*(r)^(r)A,(r, a)]] - 

N 

- Zf^T, a)Y S (r, a){T rS (r, a) + ]T 5 3 (a - ^(r))^ - e^*(r)^(r)^(r, a)]} » 0, (119) 

i=i 
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where 

T„{t,3) = -h-L**g fgl[ > - y^Y^F^F s ,)(r,a), 
2 2 

T tS (t, a) = F s{ (t, 5)7r f (r, a) = [B(r, a) x tt(t, a)],, (120) 

are the energy density and the Poynting vector respectively. 

Since the canonical Hamiltonian is (we assume boundary conditions for the electromag- 
netic potential such that all the surface terms can be neglected; see Ref. 0) 

N r 

H c =-J2 K i?( T )vK T ) + I dV[7r A (r, o)d T A A (r, a) - PfM (r, a)z?(r, a) - C(t, B)\ = 

1 = 1 

= J d 3 a[d f (-K f (r, a)A T (r, a) - A T (r, a)T(r, £)] = -/ d 3 aA T (r, a)T(r, a), (121) 

with 

jv 

r(r,a) = -y/fr.^+^e^M^M^^-^r)), (122) 

i=i 

we have the Dirac Hamiltonian (A At (r, ex) and A T (r, a) are Dirac's multipliers) 

F = J d 3 a\»(T, a)H^T, a) + A r (r, a)7r T (r, a) - A t (t, a)T(r, a)]. (123) 

The Lorentz scalar constraint 7r r (r, c?) is generated by the gauge invariance of S; its 
time constancy will produce the only secondary constraint (Gauss law) 

r(T,?)»0. (124) 

The six constraints Tl^T, a) ~ 0, 7r T (r, a) ~ 0, T(r, a) « are first class with the only 
non vanishing Poisson brackets 

{ft M (r,a) , W v (t,ct)} = 

= {[/^(r, a)z fv (r, a) - l u (r, a)z ffl (T, a)\ . h 

V7(t, of) 

+ ^(r, <?)F„(t, a)^(r, a)}T(r, a)5 3 (a - a') « 0. (125) 

Let us remark that the simplicity of Eqs.( |125] ) is due to the use of Cartesian coordinates: 
if we had used the constraints TLi(t, a) = P(r, (r)T-C^(r, a), Ttf(T, a) = z?(t, <t)7-^(t, a) (i.e. 
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nonholonomic coordinates), so that their associated Dirac multipliers A/(r, a), A r -(r, a) would 
have been the lapse and shift functions of general relativity, one would have obtained the 
universal algebra of Ref . [|TJ . 

The ten conserved Poincare generators are 

P» = J d 3 aff(r,a), 

Jr = / d 3 a (z"(t, a)p»(r, a) - z»(r, a)ff(r, a)), (126) 

so that the total momentum is built starting from the existing energy momentum densities 
on the hypersurface 

I d?aW{r, B)=tf- j d 3 al„(r, a)[T TT (r, a) + 

N 

4 = 1 

Vi yJrn%-Y s {r, a)[K if {r) - e^(r)^(r)A,(r, a)]Mr) - e^*(r)^(r)A,(r, a)]] - 
- J d 3 az fl ,(T, a)Y S (r, v){T tS (t, a) + 

N 

+ Y, - Vi(r))[K is (r) - e^(r)^(r)A,(r, a)}} « 0. (127) 
i=i 

If we add the gauge- fixings (flTf ) and the Dirac brackets (|5T|), fl57"D, we remain with the 
variables x^,p^, b^ Sg u , A^, n A , ffi, Hi, 9i, 9* and the twelve constraints 

HUr) = J d 3 aW(r, a) = pP - J d 3 a[7r 2 (r, a) + B 2 (r, a)] + 

N I : 2 

+ E^V^ 2 + - e i eUr)e i (r)A(r,ff l (r))] } - 

i=i 

r - N 
- b^(r){ / d 3 a[B(r, a) x tt(t, B% + CO - e^*(r)^(r)A f (r, ^(r))]} w 

J i=l 

H^(t) = #( r ) J d 3 aa f H u (r, a) - b u f {r) J d 3 aa f H»{t, a) = 

= Sr(r) - mr)K - b»(r)K] [1 J d?aa f [vf 2 (r, a) + B 2 (t, a)] + 

N l = 

+ EtfOW™? + [*O0 " eMr)ei{r)A{rMr))] 2 ] ~ 
i=i 

~ mmir) - ^(r)6j?(r)] [ / cfW [J3(r, a) x 7r(r, a)] s + 
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N 

i=i 

vr T (r, a) « 

r(r,a)^0, (128) 

with Poisson algebra 

{W{r),ir{T)Y = J <Pa{[Wf{r) - W(t)]tt,(t, a) + 
+b?(T)F f . s (T,m(r)}T(T,a) 

{W"(T), W^(r)}* = / dW{[&?&?(r) - 6£6?(T)]tf(T)7T,(T,*) + 

+6^(r)F M (r,a)[fo?(r)6f(T) - 6?(r)6?(r)]}r(r, a) 
{H^{r) } n aP {r)Y = C^ aP W 5 (r) + J d 3 aa*a* 
{6£(t)&?(t)([6^(t) - b?Mr)]Mr, 3) + 6?(r)F„(T, *)6?(t)) - 
-6g(r)6f (r)([^6?(r) - W(T)]7r, (r, a) + 6£(t)F„(t, *)6?(t)) - 
-6£(T)6J(T)([6Jf6?(r) - 6^(r)]7r,(r,a) + 6jf(r)F„(r, *)&f (r)) + 
+6SW6?(r)([W(r) - W(t)Mt, a) + #(t)F„(t, (r))}r(r, a), (129) 

and the form of Eq.([5(J) for the Poincare generators. 

Then we make the canonical transformation of Eqs.([59"D, so that the new variables are 

x%,p%, Sp, A A , vr A , r/ 4 , « i5 U 9* and one has 

S? B = e A ^u{p s ))e B v {u{p s ))S^ « [#(r)&? - 6f (r)tf] [± / [tt 2 (t, a) + 5 2 (r, a)} + 

N 



i=i 

+ P#(r)6f (r) - & r B (r)^(r)] [ / dW [B(r, a) x tt(t, a)], + 



N 

+ T,Vl(rM(r) - e^*(T)^(r)A s '(T,^(r))] ]. (130) 
i=i 

The final gauge-fixings d6^) reduce the variables to x^,p%, A T , 7T T , A r , 7T r , /yj, fvj r , 6?j, 9* , 
with "r" being a Wigner spin-1 index and o = r a Lorentz-scalar one. 
We have 
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* / = 

i=i 

+ fafr? - %V) [ / dW [5(r, if) x 7f(r, a)], + 

N 

+ E^)K%) - eA*(r)^(r)^(r,^(r))] ] 

8=1 
2V 

« E(^)kV) - e^;(r)^(r)^(r,^(r))] - 
i=i 

N 

- E^MKM - ^4^))] ) - 

1=1 

- y dV (<r r [B(r, a) x tt(t, a)]" - cx s [5(r, a) x 7?(r, a)]"; 



iv 



Sf » -Sf = -E^HV™? + fe(r) - e i e ( *(r)e < (r)A(r,^(r))] a - 
i=i 



./;" ^ ■ ■ tA=. (131) 

Ps+Vs^/Pj 

The Poincare generators now have the form of Eqs.([7TD and only six first class constraints 
are left 



H»{t) = - t*"(u(p.)) [~ J d 3 a[n\r, a) + 5 2 (r, <?)] + 



+ E W™f + - eA*(r)^(^)^(r,^(r))] 2 ] - 



- e{f(u(p.)) [- / d 3 a[B(r,a) x tt(t, a)] + 



N 

+ EK( r ) " e^*(r)^(r)^(r, #(r))] ] » 
i=i 

7r r (r,a) « 



r(r, 5)«0, 



= / dV{[^(«(p s ))e^( U (p.)) -^ T (u(p s ))e^u(p s ))}ir r (r,a) + 
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+e?(u(p s ))F rs (T, ay s (u( Ps ))}T(T, a) (132) 



or 



N 



W(t) = Vs^Jpl ~ + - e^I(r)^(r)A(r,^(r))]2 + 

+ ~ y rfV[vf 2 (r, a) + 5 2 (r, a)] « 

AT 
i=l 

- /" dVB(r, a) x 7? (r, a) « 
7r T (r, ct) « 

r(r,5)w0, (133) 

Then, for N > 2, we do the canonical transformation (|73"D, ([741): the new form of the 
constraints is 

N 

n(r) = e s - Q> X 

i=i 



I AT-l ^ N-l 

\ ™< + 1^k + (t) + *fA(j) - 6^9*^)9^) A(r,ff + (r) + -= ]T w (r))P 

N iV a=l ViV a=1 

+ i y rf 3 a[vf 2 (r, a) + £ 2 (r, a)]} = e, - S (P+/)s - E (F)s « 



JV ^ AT-l 



W p (r) = £ + (r) - E^(r)^(r)A(r,i7 + (r) + ^ E 7a*Pa(r)) 



i=l V 1 v a=l 



- y d 3 aB(r, a) x 7? (r, a) = P(p+i) s + P(f)s ~ 
7r T (r, a) « 

r(T,a)wO, (134) 

where -E^s = |/ <i 3 cr[7r 2 (r, a) + -B 2 (r, a)] and -P(f)s = ~I d 3 <J7i(T,a) x B(r,a) are the 
rest-frame field energy and three-momentum respectively [now we have 7?(r, a) = -E(r, a)], 
while E(p + i) s and P(p+/) s denote the particle+interaction total rest-frame energy and three- 
momentum, before the decoupling from the electromagnetic gauge degrees of freedom. 
The final form of the rest-frame spin tensor is 
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JV i JV-l 

S7 = -Y,^Kr)9 l (T)[ V l(r)A s (r,ff + (r) + -= £ %p b (r)) - 



i=i v -/V 6=1 

i JV-l 

- n s + (r)A r (r,ff + (r) + -= ^ %p b (r))) + 

JV-l i JV i JV-l 

+ E K(r)«(r) -tsEwWW^+W + "TXF E 76iP6(r)) 
o=i v iV i=l ViV 6=1 

1 JV i JV-l 

- P 8 a (r)K(r) - - 7W Y.1a^ d t{T)e i {T)Ar{T 1 ri + {T) + -= E 7hA(r))] - 

VN i=1 ViV 6=1 



J d 3 a (a r [B(t, a) x n(r, a)}* - a s [B(r, a) x tt (r, <?)] r ) = + 5, 



rs 

(F)s 



N 1 JV-l 



i=l viv a=1 



i JV-l j JV-l 

. m? + [^k+(t) + v/iV E 7 ai v? a (r) - e^*(r)^(r)l(r, r/ + (r) + -= £ 7«Pa(r))] 2 

-~Jd 3 a<T r [n 2 (T,a) + B 2 (T,a)}, (135) 

while the Dirac Hamiltonian is 

H D = \{t)H - \(t)H p + J d 3 a[X T (r, a)n T (r, a) - A t (t, o)T(t, a)]. (136) 

We can check that 0*(t)0j(t) is a constant of motion for each i=l,..,N, so that we will write 
it as 9*9i. 

Let us now look at the electromagnetic Dirac observables, namely at the functions on 
phase space invariant under electromagnetic gauge transformations. Referring to Ref. for 
the detailed calculations, for the electromagnetic field we have the following decompositions 

d 

A r (r, a) = g^Ve m (r, a) + A r ± {r, a) 

1 d N 
tt p (t, a) = vrKr, a) + £^[T(t, a) - E ^^{9 - ffi(r))) 

Vem(r,a) = -±-^.A(r,a) 



{ Vem (r,a),T(r,a')r = -S 3 (a-a') 
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K ± (r, a), vrKr, a')}** = -(S rs + ^)5\a - a'), (137) 

with the pairs of conjugate variables A t (t, a), ti t (t, a) ~ 0, i] em (r, a), T(r, a) spanning 
the Lorentz-scalar gauge subspace of phase space, and with A±(t, a), n±(r, a) [d a -A±(T, a) = 
d a ■ t?±(t, (?) = 0] being a canonical basis of electromagnetic Dirac observables, transforming 
as Wigner spin-1 3- vectors. 
Since we have 

K(r),r(r, <?)}** = 

K(r), T(r, a)}** = -e&Bi^ia - i*(t)), (138) 

the particle momenta /%(r) are not gauge invariant. Also the Grassmann variables 
0j(r), 0*(t) are not gauge invariant 

{^(r),r(r )C ?)r = -^(r)5 3 (a-^(r)) 

{^(r),T(r,a)r =ie i 91(T)5 3 (<?-fii(T)). (139) 
The Grassmann Dirac observables are 

0.( r ) = e-^^W)^^), {^( r ), r(r, a)}** = 0, 
0*( r ) = e ^(^M)^( r ) ; (0*( r ), r(r, a)}** = 0, 

0$ = 0*0,, r(r, a)}** = o. (140) 

Since we have 

iu(r) -e^eJir,^)) = ^(r) - e^A^r, j£(r)), (141) 

with 



K(r),ft5(r)r = «5y6", 

K(r),^(r)r = 0, K(r),fl;(r)r = 0, (142) 
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the particle Dirac observables are ^(t), /%(r): the scalar particles have been dressed with 
the Coulomb cloud (like the fermion fields in Ref. 
From Equations (|137|) , we have 

Id N - - 

n\r,a) = {jf ± ( T ,3) + — — \T(r,ff) - £ e^ 3 (<? - ^(r))]} 2 « 

1 d 

« 7?i(r, a) - 27?j.(r, <?)•——£ e, ^(a - jft(r)) + 

A CT 9a ^ 

TV 1 Q 1 Q 

+ £ M^]M;^a;^ 3 ^-^W) • a:^( ? -^W). ( 143 ) 

so that by integrating by parts and using the property Qf — 0, Qi — ei9*9i, of Grassmann 
variables, we obtain 

J d 3 an 2 (T, a) = 

1..N x 

= j d 3 a{nl(r, a) - £ Q % Q 3 5\a - ^ T ))—5\a - ff 3 (r))} = 

. 1..N 

= j dVLfo a) - £ c(^(t) - ^(r)), (144) 
where we have introduced the Green function c(x) defined by 

= 4fi£rr (145) 

Having decoupled the electromagnetic gauge variables, the canonical basis of Dirac ob- 
servables x%(t),p%, A_l(t, a), 7t_l(t, a), ffi(r), Ki(r), 9i(r), 9*(r) spans the phase space, where 
the remaining gauge invariant four first class constraints in Eqs.( |134] ) have the form 

n(r) = e. - {£th 

i=l 



-y N-l y N-l 

x ,af + [^k+{r) + v / iV£ %A{r) - e^M-ifa ff + {r) + -= £ %ip a {r)W 

\ iV a=l ViV a=1 

i 1..N i JV-1 

1 i+j viv a=1 

+ \j dM#L(r, <?) + 5 2 [A ± (r, a)] ]} « 
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N 1 N-l 



n p (r) = k+(t) - ^e^*M±(r,r/+(r) + -= E %iP a {r)) - 
»=i V iv a=1 

- y d 3 aB[4(r, 5)] x t? ± (t, a) « 

{ft,7^}** = {7^, 7^}** = 0; (146) 

the Bianchi identity d ■ B(r,a) = has been used to obtain Hp. The variables K+,7c a are 
the Dirac observables induced by «j. Now the Dirac Hamiltonian is H D = A(r)7Y(r) — A(r) ■ 

The rest-frame spin tensor becomes 

N i JV-l 

i=i viV b=1 

i AT-1 

- V Ur)A r ± (T,ff + (T) + -= + 

V iV b=1 

N-l i AT i AT-1 

+ E [PaM(^) " -7=Ewfi?Mi (r,ff + (r) + -= E 7taA(r)) - 
o=i viv i=1 ViV b=1 

i AT i AT-1 

- P s a (r)(K(r) ~ -7=Ei.4M r i(M;(r) + -7= E 7^))] - 

ViV i=1 ViV b=1 

N ~ ~ r Id 
-^0*0, J d 3 a(* r e suv -a s e™)[— ^-5 3 (a - ^(r))} B v [A ± (r,a)} - 

- J d 3 a (a r [B[A ± (r, a)} x tt ± (t, a) f - <r s [S[A ± (r, a)] x tt ± (t, a) f) = 

N-l 
a=l 

- I (<7 r [S[l ± (r, a)] x tt ± (t, a) ]' - a s [B[A ± (r, a)} x tt ± (t, a) f) - 

N 1 A"-l 

-E^mV r + (r)A s ± (r,ff + (r) + -= E TtaAOO) - 
i=i viV b=1 

1 AT-1 

- ^(r)A;(r,77 + (r) + -= E TbiAOO)] + 

viV b=1 

j AT-1 AT j ^-1 

" ^7 E E Taiei^MKr, V + (t) + -= E 76iP 6 ( 

ViV a=1 i=1 ViV b=1 

N j AT-1 

- Pa(^) E 7aM*Ml(r, 77+(r) + -7= E 7WP&0"))] - 

i=i viV b=1 



r)) 
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e^^em / d 3 ac(a-rfi(r)){2B k [A ± (r,a)] - 
i=i J 



da 

°s,FREE °s,INT 



Cjor cjr 



N 



N-l 



%iPa( T ))Vi X 



N-l 



N-l 



o=l ViV a=1 



1 r 1 <9 ^ 

- ] rfW [-27f ± (r, a) ■ — — £ em.5\a ~ *(r)) + 



A? 



' ^W[7fi(r,5 ! )+ J B 2 [A ± (r,5 ! )] 



N 



N-l 



i=l 



\ 



j iV-1 j JV-1 

mf + [-^+(r) + v/iVE Tai#a(r) - e^*M±(r, r/ + (r) + -= ]T 7«p a (r))] 2 



a=l V J » a=l 

JV L.jV i a 

+ E e M(E ^^^[7— TTT^Kr) - %(r)) - I7j(r)c(^(r) - %(r))] + 



8=1 



(147) 



Therefore the final result is the extraction of the static action-at-a-distance Coulomb 
potential from the electromagnetic field theory, which is reduced to the transverse radiation 
field. This is obtained at the pseudoclassical level, where the hypothesis of charge quanti- 
zation is reflected in the property Q\ = 0, which regularizes the classical electromagnetic 
self-energy and produces the rule J2i^j without the need to introduce Feuynman- Wheeler 
54fl theory of the absorbers. Clearly all the effects of order ef (like those producing the run- 



away solutions in the Abraham-Lorentz-Dirac equation; see Refs. |55|J56[| ) are killed: only 
the interference effects of order e^ej, i 7^ j are preserved. These problems will be investi- 
gated in more detail elsewhere. In any case, it is clear that in this formulation the initial 



63 



data problem consists in specifying Cauchy data for the radiation field and Newton-like 
action-at-a-distance initial data for the charged particles with mutual Coulomb interaction. 

In particular, if we put equal to zero the radiation field, A±(r, a) = n±(r,a) = 0, we 
obtain a well defined relativistic 1-time reduced N-body problem with the Coulomb potential 
extracted from field theory and with a decoupling of the motion of the center of mass, whose 
constraints are 



N 



H(t) 



e s - {E^\ 
i=i \ 



N-l 



K+(r) + VN ^ %Mr)] 2 ~ 



a=l 



1 UN . . I 

2 E MttlM^i] c[-y= £ (%i - l aj )Pa(r)]} 



N-l 



a=l 



N 



{E^\ 
i=i \ 



N-l 



m: 



a=l 



i 1..JV i N-l 

o E Mi^i][cM] c[-= E (7- - 7 aj )Pa(r)]} « 



W p (r) = k+(r) ^0 



{w, w;}** = {% n;y* = o ; (148) 

whose natural gauge-fixing for decoupling the last three constraints is ff + (r) ~ as in the 



free case [the nonrelativistic limit, using the methods of Refs. | lq| , would give the first class 
constraint Ti, = Er — Hr pa (where Er and Hr are the reduced energy and Hamiltonian 
respectively) for the parametrized Newtonian N-body problem with Coulomb interaction]. 



See also Ref. [|51|] and its bibliography for models with additive potentials like the Coulomb 



one in Eq.( |148[ ). For the N=2 case, Eqs.( |148| ) were already found in Refs. |57 |, where there 



is a study of the classical equations of motion and their explicit solution in the equal mass 
case. 

It is difficult to find the gauge-fixings replacing ff + pa 0, i.e. which is the point of the 
system to be identified with x%, for the full theory, and then the final form of the rest-frame 
spin tensor (|147 ). Given these gauge-fixings, we have to find a canonical basis of Dirac's 
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observables with vanishing Poisson brackets with the constraints 7i p ~ of Eqs.( |146| ) and 
with the gauge-fixings, and to reexpress Ti. ~ of Eqs. ( |146| ) in terms of them: this will 
give the final 1-time reduced Hamiltonian for N charged scalar particles interacting through 
the Coulomb potential and coupled to the radiation field. However, an implicit form of the 
three gauge-fixings may be found by remembering that in any instant form of the dynamics 
there are only four generators of the Poincare group which depend on the interaction. In 
the rest-frame instant form they are e s = t] s ^Jj^, given by the first of Eqs.( |146| ), and the 
three boosts J°\ given in the last of Eqs.( |13i|) [they depend on the interaction through the 
term r] s J^[. The other six, i.e. p s and J' l J [namely S r s % must be identical to the analogous 
generators in absence of interactions, namely to the ganerators for N neutral scalar particles 
plus the electromagnetic field. Therefore, from Eqs. ( |147| ) we must require the vanishing of 
the interacting part S r s s INT of the rest-frame three-spin 

i N -, N-l 

X\r) = ~e krs S™ INT = -e fc ^;(r)£e^Mi(T,r7 + (r) + _= £ 7 6i p 6 (r)) - 
z i=i viv 6=1 



e krs N-l N y N ~ X 

= E P r a (r)J2"faie i 9*e t A s L (T,ff + (r) + -= £ %p b (r)) 



a=l i=l VJ« 6=1 

N . W 

-E e A*^ J d 3 ac(a-rf l (r)){2B k [A ± (r,a)] - 



i=i 



- a . —B k [A ± (r,a)] + a ■ ^B[A ± (r,a)]} « 0. (149) 

The time constancy of these constrains determines the Dirac multipliers A(r), so that 
the final Dirac Hamiltonian becomes Ho = A(t)7Y(t). The problem of evaluating the Dirac 
brackets of the three pairs of second class constraints 7i p ~ and x r ~ as well as the 
attempt of diagonalizing these brackets will be studied elsewhere. 

Let us remark that, if in Eqs. (|146|) we disregard the constraints 7i p ~ and eliminate 
the field energy term in ii. « 0, we get a constraint 7i' ~ in which Aj_ is an external 
radiation field coupled to the system. 

If we eliminate the particles and we remain with only the field, the four constraints are 

H{t) = ts~\j dV[7?i(r, a) + B 2 [A ± (t, *)]]}« 
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H p (t) = / d 3 a B[A ± (t, a)} x tt ± (t, a) = P (F)s « 



{n, n;}** = {n r p , n;}** = o. (150) 

Now the gauge-fixing should be the vanishing of a center-of-mass field three-coordinate 
conjugate to the total three-momentum P(f)s = 7~t p ~ 0. Therefore, we have indirectly 
shown that there must exist a center-of-mass decomposition also for classical gauge field 
theory and that this is the lacking ingredient to get a reformulation of classical field theory 
as a "1-time rest-frame field theory" with Wigner covariance of the variables relative to the 
center of mass when restricted to configurations belonging to timelike representations of the 
Poincare group. Again, only three center-of-mass coordinates would not be covariant. 
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VII. THE 1-TIME THEORY VERSUS THE FESHBACH-VILLARS 

HAMILTON! AN THEORY 



In Ref. J571 (see also Ref. the first quantized one-particle second order Klein-Gordon 
equation coupled to an external electromagnetic field, (D 2 + m 2 )0(x) = 0, = <9 M — 
ieA^x), was put in Hamiltonian form after having rewritten it in the first order formalism: 
<Po( x ) ~ i; D o<P{x) = 0, (D 2 - m 2 )4>(x) + imD <j) (x) = 0. If we put <f>(x) = ^[<f(x) + x(x)}, 
4>o{x) = ^[<f(x) ~ x( x )], w e obtain 

1 2 

id Q ip(x) = - — (-id - eA(x)) Up(x) + x{x)) + (eA (x) + m)<p(x) 
2m 

1 2 

id x(x) = - — (-id - eA(x)) (tp(x) + x(x)) + (eA (x) - m)x{x). (151) 
2m 

In a 2 x 2 matrix formalism we have (t* are the Pauli matrices) 



idn^ix] = H ^(x) 



y(x) 



I <p(x) ' 
\X(x) j 

1 2 

if = - — (— i5 - eA(x)) (r 3 + zr 2 ) + mr 3 + eA Q (x)l 
2m 

->A„->o H = — (r 3 + ir 2 ) + mr 3 . 

2m 



(152) 



-*2 

In the momentum representation we have H a = ^(t 3 + 2t 2 ) +mr 3 and this Hamiltonian 
can be diagonalized (p° = +^/m 2 + p 2 ) 



H , u = U- 1 (p)H U(p)=p°T 3 

^u(p) = U-\p)^(p) 
id ^u = H ^ u 



( V m 2 + p 2 ^ 







■y/m' 2 + 



i 



2^/mp 
1 



= [(m + - (m -p°)n] 



2y/mp 



= [(m + p°)l + (m-p°)n]. 



(153) 
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Like in the case of the Foldy-Wouthuysen transformation for particles of spin 1/2, also in 
the spin case the exact diagonalization of the Hamiltonian cannot be achieved in presence 
of an arbitrary external electromagnetic field |37| . 



At the classical level, the first class constraint (p — eA(x)) 2 — m 2 may be resolved 



as (p° — eA°{x) — \j m 2 + {p — eA{x)) )(p° — eA°(x) + \j m 2 + {p — eA(x)) ) ~ 0, which can 



be replaced by the 2x2 matrix 

/ I ~ 2 \ 

' (p° - eA°(x) - \Jm 2 + (p- eA(x)) ) * 



y (p° - eA°(x) + \lm 2 + (p - eA(x)f) j 



eA°(x))l + \jm 2 + (p- eA(x)) 2 T 3 w 0. 

This shows that the 2x2 matrix formalism has a topological origin in the two disjoint 
branches of the timelike orbits. However, in presence of the external electromagnetic field 
p M is no more a constant of the motion and it is not clear how to apply the theory of the 
canonical realizations of the Poincare group. 

In contrast, in the 1-time theory for N=l in presence of a dynamical (not external) 
electromagnetic field, from Eqs. (|133|) , (|134|) [in the canonical transformation (|73|) , (ff4]), ff(r) 
and k(t) remain fixed], we have the following rest-frame first class constraints 



H{t) = e s -r 1 \jm 2 + [k(t) - e6*6A(r, r/(r))] 2 --J d 3 a[7r 2 (r, a) + B 2 {r, a)} w 
H ( T ) = k(t) - e6*6A(r, ff(r)) + f d 3 an(r, a) x B(t, a) w 



7r T (r, a) ps 



T(r, a) = -d- tT(t, a) + e6*65 3 (a - ff{r)) w 0, (154) 
or in terms of Dirac observables [see Eq. (|146 ) ; again we have (ee*9) 2 = 0] 



H(t) = e s - V \lm 2 + [k(r) - eHA ± {r, ff(r))} 2 -\jd 3 a (vr^r, a) + B 2 [A ± {r, a)}) « 
7? p (r) = k(r) - e§*9A ± (T, ff(r)) + J d 3 a7r ± (r, a) x B[A ± (r, a)} w 0. (155) 

Even if it is not yet clear which are the natural gauge-fixings to decouple the constraints 
Tip ~ [i.e. which is the point of the plane to be identified with x^\ in absence of the field, 
the gauge-fixing is ff(r) ~ 0, namely the particle position is described by x%] so to remain 
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only with one final first class constraint 7t ~ for the final Dirac observables, some remarks 
may be done. 

There is a well defined coupling of the electromagnetic field to each branch 77 = ± of 
the particle mass spectrum, independently from the other: in a sense for p 2 s > we have 
obtained a diagonalization of the overall particle +field mass spectrum. However this has 
been achieved by treating the electromagnetic field as a dynamical field and not as an 
external one. Moreover, we have used a definition of rest frame based on the conserved 
Poincare generator p^, which is associated with the isolated system particle+field and which 
does not exist when the field is external. Finally, even if we can form a 2 x 2 matrix 



constraint 



' n\ v =+ » 



\ H\ n =- j 



0, both 7^1^ contain the term in the field energy: only after 



having found the rest-frame field theory based on a decomposition of the fields in center-of- 
mass and relative variables we will have the tools to try to understand better this matrix 
formulation and its connection with the Klein-Gordon equation (the coupling of the 1-time 
theory could be a nonminimal nonlocal coupling of the original Klein-Gordon equation to 
the electromagnetic field). 



In Ref. there is also a discussion of the first-quantized effect of zitterbewegung of 
the canonical noncovariant Newton- Wigner 3-position operator x op for a scalar particle and 
of its nonlocalizability. Since the positive energy states (with the standard scalar product) 
do not satisfy a completeness relation, the narrowest wave packet containing only positive 
frequencies has a width of the order of the Compton wavelength 1/m. Equivalently, the 
eigenstates of x op = id /dp contain both positive and negative energies: in the diagonal 2x2 
representation we define a "mean position operator" [which has 5 3 (x — x) as eigenf unctions] 

as x UjOP = U- l (p)i^U(p) = i-§~l- 2 {rJ+^) T ^- = + %i where the odd matrix r i 
couples positive and negative energies. The operator x^' , with velocity doX^) = —r^= — t 3 , 

' V m 2 +p 2 

describes the steady motion of the particle, is the Newton- Wigner position operator in this 
representation but has eigenfunctions delocalized over a region of radius of the order 1/m 
(they have the previous velocity as group velocity); instead a wavepacket with a localization 
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more accurate than 1/m (built with eigenfunctions of xu,op) will contain both kinds of 
energies and will have a trembling motion (zitterbewegung) superimposed to the steady 
motion in a region of the order 1/m. 

In Klein-Gordon quantum field theory ||59|| , the byproducts of these effects are 
the nonexistence of perfectly localized states (eigenstates of a local occupation num- 
ber operator Ny }0p (x°) = J v d 3 x N op (x, x°) for an infinitesimally small volume V; 
[Nv,op(x, x°), Nv, op (x , x°)\ is different from zero for distances | x — x\ of the order 1/m 
we cannot fix the particle number in a volume smaller than 1/m and the number of particles 
in two volumes V and V' cannot be measured simultaneously unless V and V' are separated 
by distances greater than 1/m. In contrast, we can define a local charge operator Qy and 
build states with perfect localization of charge in infinitesimally small volumes; however 
[Qv, N op ] ^ 0, i.e. the eigenstates of Qy do not have a sharp particle number. 

In contrast, in the classical rest-frame 1-time theory of a scalar particle, positive and neg- 
ative energy states are completely decoupled (the upper and lower branches), and, instead of 
the zitterbewegung effect, we have the impossibility to localize the center-of-mass canonical 
3-position x s [to which collapses the particle position after the gauge-fixings ff(r) ~ 0; the 
first class constraint becomes 7i = e s — 77m « and we have a Hamilton- Jacobi description 
of the free particle, with z s /r]m as Jacobi data] inside the noncovariance worldtube of radius 
| S | /m in a frame independent way. 
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VIII. 1-TIME RELATIVISTIC STATISTICAL MECHANICS 



In Ref. |5T |) there is a a review with rich bibliography till 1966 of the two types of 
approaches to classical relativistic statistical mechanics: i) the nonmanifestly covariant and 
ii) the covariant ones. While the former approaches started in 1911 with Refs. ||62|| , where 
the equilibrium of a relativistic (classical Boltzmann and also quantum either Bose or Fermi) 
gas was studied, only in 1940 in Ref. a first relativistic Boltzmann equation was given 
in the framework of relativistic kinetic theory; then the studies in plasma physics focused 
on equal-time Hamiltonians 

N I Z 2 1 r 

H = Y,{\jm 2 l + (p i -e l A(x l ,x°)) +V(x h x°)} + - / d 3 z[n 2 (z, x°) + B 2 (z, x°)], (156) 

i=i 4 J 

derived a relativistic Liouville equation for a density involving both particles and fields and 
were extensively studied in the sixties [see the bibliography of Ref. |6l| . 

In contrast, the covariant approach started with Ref. [Q, where, taking into account 
constraints, there are the first steps in defining covariant relativistic kinetic theory, relativis- 
tic statistical mechanics (canonical ensembles) and relativistic thermodynamics following 
the first developments of relativistic hydrodynamics [see the bibliography of Ref. fl6l|1; the 
Vlasov-, Boltzmann-, Landau-, Fokker- Planck- relativistic kinetic equations came out as ad 
hoc semiphenomenological equations [again see the bibliography of Ref. J5T]]]. 

In Ref. ||61|| ) there is a critical analysis of the problems in the foundations of relativistic 
kinetic theory and relativistic statistical mechanics stemming from the covariant relativistic 
description of particles with either action-at-a-distance or field-mediated interactions. Due 
to the absence of an absolute time, relativistic kinetic theory (//-space) is a statistical theory 
of curves (worldlines) and not of points (Cauchy data) as at the nonrelativistic level. Anal- 
ogously, due to the presence of N times in the description of N particle systems, relativistic 
statistical mechanics (r-space, Gibbs ensembles) was a statistical theory of N-dimensional 
manifolds (the gauge orbits spanned by the N times in the 7N constraint manifold defined 
in the 8 N-dimensional phase space by the N mass-shell first class constraints p 2 — m 2 « 0) 
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instead of 1-dimensional manifolds as in the nonrelativistic case (so that one has N and 
not 1 relativistic Liouville equations in T-space and associated problems with the defini- 
tion of measures). These problems were already posed in Ref. ||64|| . However, there was a 
not clear understanding that the transition from the configuration space [worldlines with 
coordinates qf (r) or qi(q°); see predictive mechanics |J5| and the predictive conditions on 



the relativistic forces for the implementation of a realization of the Poincare group |66]] to 



phase space [with coordinates xf(r),pf(r)] had to face the No-Interaction-Theorem []13| , [T4 
so that gf(r) = xf (r) is allowed only in the free case. Moreover, in Ref. [61] there is the 



recognition that to build relativistic kinetic theory and statistical mechanics the initial data 
must be chosen as the mathematical Cauchy data on a spacelike hypersurface [it is not clear 
whether they can be measured] and not for instance as data on the backward light-cone of an 
observer [in the case of the electromagnetic field, the light cone is a characteristic surface of 
Maxwell equations and does not define a well posed Cauchy problem: the future field is not 
completely and uniquely determined]: in this case, in /x-space, starting from the predictive 
equations of motion in the first-order formalism [q^(r) = w m (t), m-u M (r) = F M ], from the 
density w M , q%, w^, r) = <5 4 [g M — g M (r, q a , u )] x 5 A [u^— w m (t, q Q , u )] and from an ensemble 
of such systems with a random distribution D (q ,u ) of Cauchy data, it is possible to de- 
fine the density in /x-space -D(g M , u^, t) = J* -R(g^, g^, w^, r)D (q , u )d 4 q d 4: u 0) to find the 
one-particle relativistic equation for D [and then, by adding ad hoc collision terms, various 
kinetic equations], to show that Af(q^,u^) = f*™ drD(q^,u^,T) is a distribution function 
such that j^{q) = J d i pM{q^ l . l p^ = mu' 1 ) u^2m9{p°)5 i {p 2 — m?) [if one chooses u 2 = 1] is the 
particle four-current, that J\f(q fl ,u fl ) satisfies the one-particle Liouville equation, which in 
turn, by adding ad hoc collision terms C(Af), becomes a kinetic transport equation. When 
the external force vanishes, the equation C{M) = is satisfied by the relativistic version 



of the Maxwell-Boltzmann distribution function, i.e. the classical Jiittner -Synge |67 
distribution [K2(x) is the modified Bessel function of order two] 

W.P^) = a Wt m 6 "^- ( 157 ) 
4:Txm 2 K 2 {mf3) 
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In this equation (3^ is the reciprocal temperature timelike fourvector |>7 \J3 = y/JP =1/ kT 
with k being the Boltzmann constant]. 



In Refs. ||68|| , there was a further development of relativistic kinetic theory and a dis- 
cussion of how to formulate covariant relativistic thermodynamics [see Ref. [^] for the 
bibliography on the subject and for a review of the associated problems], whose final for- 
mulation [containing the solution of the problem of having causal propagation of heat flow 
(parabolic transport equations replaced by hyperbolic ones) with some molecular speed less 
than the velocity of light] was given in Refs. |70| for fluids both in equilibrium and out of it. 



The conclusion was that at equilibrium all relevant fourvectors [the reciprocal temperature 
P^, the entropy four-current s M , the volume fourvector V 1 , the number density of particles 
n^, the four-momentum p M ] should be all proportional, through their rest-frame values, to 
the hydrodynamical four- velocity [u 2 = 1]: (3^ = (3u^ = u^/kT, = Vu^, s M = su^, 
= nil 1 *, p^ = w M is a timelike eigenvector of the energy- momentum tensor T^ u if it 



satisfies certain positivity conditions f67[j . 



At equilibrium, starting from the Jiittner-Synge distribution, we can derive the co- 
variant canonical partition function Qat(V^, /3 m ) for an ideal classical relativistic Boltz- 
mann gas of N scalar particles of mass m (and going to the quantum level of Bose and 
Fermi ideal quantum gases) and then the covariant grand-canonical partition function 
S(V^,/? M ,^) = J2n L Q V(V^ ', /? M ) [i M = with /i the relativistic chemical potential, 
related to the nonrelativistic one by fi = yUjvij+m] and the covariant thermodynamical poten- 
tial f2(^/ M , i^) = —kTln'Eiyt 1 ,^,^). This has been done in Ref. ffl] following previous 



results in Refs. |72],73|, after the identification of the invariant microcanonical density of 
states 



N N 



Mp 2 ,V 2 ,p,Vn = J ^(p-E^n^CPi.m), (158) 



i=l i=l 



based on the invariant momentum space measure [ff^,f75 



dcr(p, m) = 2V l ,p' 1 6(p )5(p 2 - m 2 )(Tp. (159) 
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Then the canonical partition function for the ideal Boltzmann gas of N free scalar par- 
ticles is [1/iV! is the Boltzmann counting factor for identical particles] 

Q N (V^^) = ±j d 4 pe-^a N (p\V\p,V») = 

' d*p6{p°)5{p 2 - m 2 )p i y i e- p " p »] N = 



N\ l (2tv) 3 

N\l\ [ (2ir)3J Vl J 
1 Vm 2 

= M [ 2^ K2{m(3)]N > P = 1 / I * r > ( 16 °) 
where the last evaluation has been made in the rest frame of the volume V. 

This same result may be obtained |7B| from the phase space description of N scalar 
particles with N first class constraints, by using a Faddeev-Popov measure in which N gauge- 
fixings [of the kind pi ■ Xi/yp^ — r< ~ 0] have been added to the N first class constraints 
<pi « to eliminate particles' times and energies (the many-time problem). The resulting 
form of Qn for a system of N interacting particles is [p M = Y^iLiPi] 

Qn(VW) = — J yti e-P*> l[e(p°)5(<t>i)S(Xi)det \ {^ Xj }\ (161) 

where 4>i = p 2 — rn 2 — . . . « [the dots are for the (generically unknown) interaction terms in 
the N-time theory] are the original first class constraints and Xi ~ are the N gauge-fixings. 
In Ref. [TfJ, there is also a study of a gas of N 2-body bound states based on the model of 



Refs. 44.45 



Moreover, Feynman nonrelativistic result JF7]] that the partition function of quantum 
statistical mechanics can be obtained by evaluating the density matrix for an N particle 
system from the continuation to imaginary times (x° — > —%h j kT) of the path integral used 
for the evaluation of the kernel of the evolution operator, may be extended to free relativistic 
scalar particles (p 2 — m? w 0) in the proper time gauge and the result is again J7S|. 



See for instance Refs. Jnj for applications of relativistic statistical mechanics to rela- 
tivistic plasmas. 

Coming back, after this sketchy review, to the 1-time description of N relativistic scalar 
particles, we see that now we have a natural kinematical framework (when the system is 
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isolated and the total momentum is timelike) for developing an 1-time covariant formulation 
of relativistic statistical mechanics of systems of interacting scalar particles along the lines 
of the nonrelativistic one; it could be termed "rest-frame instant form relativistic statistical 
mechanics" . 

To arrive at this formulation, let us start from the ideal gas in the formulation of Eq. ( |161[ ) 



and let us try to reproduce the canonical partition function of Eq. ( |160| ) with the new formu- 
lation. Let us make the canonical transformation from the basis xf,pf of the N-time theory 
to the basis x ,p M , T RcL} e Ra , p a , if a of Eqs.(|3ll) and then to the basis T, e, z, k, T Ra , e Ra , p a , ^ a 
of Eqs.(j35|). We have 



N N-l 



Yl d A Xid A pi = d A xd 4 p Y[ dT Ra de Ra , d 3 p a d 3 iT a = 

1 = 1 (2=1 

N-l 

ded 3 kdfd 3 z J[ df Ra de Ra d 3 p a d 3 ii ai (162) 



a=l 



since d 4 x = \jl + k 2 dfd 3 z/\ e| 3 and d A p = \ efded 3 k/\/l + k 2 ||. 



The gauge-fixings needed to go to the rest-frame instant form are T — r w and T Ra w 0. 
Since we do not know explicitely the inverse canonical transformations, we could expect 
to have a Jacobian coming out in the following transformation [to avoid the problem of 
degenaracies we start by making calculations with different masses mi\ 

N N-l 

II *(P?)*(P? - ™*)5(f(x,p) - r) J] ^f^p)} = 

i=l a=l 
N N-l 

= j n e^sip 2 - ...)s{f - T ) n 5(e Ra )5(t Ra ) = 

i=l a=l 

N N-l 

= J J] 0^)S(e 2 - ...)5{f - r) J] S(e Ra )S(f Ra ). (163) 

i=l a=l 

We have written p 2 — ... = e 2 — .. ~ for the (explicitely unknown) mass spectrum constraint; 
but with the help of Eq.flS2|) and due to the 9(rji), Eq.( |163|) may be rewritten as 



TV 

i=l \ 



N-l z N-l 

m 2 + N(Y / larfa) 5(f - r) n SifnaMen*), (164) 

a=l a=l 

where we have put rrii = m, because now there is no ambiguity with the uppest positive 
branch of the mass spectrum. 
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By using Equations (|162|) , ( |164j ), and our gauge- fixings, the analogue of Eq. (|161|) is [the 
Jacobian J is to be fixed by comparison with Eq, ( |160| )1 

i „ N N-l 

Qn{V»,P) = —J ^ Ue-^e(p°)5(p2 - m 2 )5(f (x k , Pk ) - r) J[ 5(f m (x k ,p k )) 



N d 4 rd 4 n 
II 



\ (2vr)3 
1 J ■ 



N\ (2Tx f N Jv» ^ § \ 



N-l z 7V-1 N 

m? + N( ^2 7a,vf a ) )5(f - r) n ^(^ a )5(e Ra ) I] e"^ 

a=l a=l j=l 



JV-1 



. JV-1 N 

■ J v d 3 zd 3 kl[d 3 Pa d 3 7r a l[e-^ p i. (165) 



ded 3 kdTd 3 z dTR a deR a d 3 p a d 3 Ti a 

a=l 

^ J <• 

iV! (2tt)^ a* 

What is lacking in the N-time theory is the evaluation of , because we do not know 
explicitely the inverse canonical transformation. 

Let us shift to the 1-time description on the rest-frame hyperplane: in the rest 
frame one had the constraints k + ~ and the gauge-fixings ff + ~ [identifying the 



rest frame center of mass with x%], so that Jm 2 + R 2 ~ yfn 2 + N(J2 a= i lai^a) an d 



~ (3ym 2 + N(J2a=i lai^a) due to Eq. flgQ|) . But if we look at the rest-frame hyper- 
plane from an arbitrary frame in Minkowski spacetime, the Minkowski canonical 3-position 
and 3-momentum of the center of mass are z s /e s and e s k s respectively with d 3 z s d 3 k s = 



{d 3 z s /\ e s \ 3 ){\ e s \ 3 d 3 k s ) [it is a Lorentz scalar @: d 3 z s d 3 k s = (Jl + k 2 s d 3 z s ){d 3 k s / Jl + k 2 s )} 



and we would like to identify d 3 zd 3 k = d 3 z s d 3 k s with the Lorentz scalar one d 3 rj + d 3 K + . In 
this way ff+ and R + would be Wigner spin-1 quantities simulating on the rest-frame hyper- 
plane [relaxing the constraints k + ~ 0] the motion of the three-dimensional center of mass 
in an arbitrary Minkowski frame: the 3-coordinates and the 3-momenta of the particles in 
this descritpion would be the Wigner spin-1 3- vectors ffi and «j in Eq.([F3|). 



This suggests the replacement of the unknown /3^pf with f3ym 2 + k 2 , which is a Lorentz 
scalar [f3 is the rest-frame reciprocal temperature]. With this prescription, by using 
d 3 r} + d 3 K+ Yla=i d 3 p a d 3 7i a = YliLi d 3 r]id 3 Ki and by replacing with its rest-frame value 
V = J v d 3 r]i, Eq. (|165|) becomes 
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1 T N r , 

= — TTr V f ^3 -/3^+^1 = 



i=i 

1 J r 2V 



N 



N\ 2 7vL (2tt) 3 

= j^sS^ 1 " (166) 

Therefore, if we choose J = 2^, we get the following definition of the Lorentz scalar 
canonical partition function of an ideal Boltzmann gas in the rest-frame instant form 

19 r N 

Qn(V,(3) = ^(^) W / y Il^^ 3 ^-^. (167) 
If, following Ref. f76|, we consider a relativistic gas formed by N 2-body bound states 



defined by the following 2-body reduced Hamiltonian [obtained from Eq.( |lll|) with the 
previous prescriptions and where Ri 2 = rji — f] 2 ] 



H R = \jm\ + V X {R\ 2 ) + k\ + \jm\ + V 2 (R 2 12 ) + k 2 2 + U(Rj 2 )' (168) 
then the definition of the Lorentz-scalar canonical partition function is [ffi 2 = \(fj\ + fji)] 
Qn(V,(3) = {2n) ] Nm [ J y e-^dSxd^^^f = 

1 A [f e-^d^R^n.d 3 ^. (169) 



{2tt) 6N N\ l Jv 

Clearly this can be extended to more complicated situations. When a better understand- 
ing of the gauge-fixings for the case of charged particles interacting with the electromagnetic 
field will be obtained, Eqs. (|146|) will serve to define the covariant analogue of the canonical 



partition function associated with the noncovariant Hamiltonian of Eq. ( |156|) . Finally, Feyn- 
man path integral approach should now be based on the continuation of the scalar rest-frame 
time T s to imaginary values —ih/kT. 
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IX. CONCLUSIONS 



In this paper a covariant 1-time formulation, the rest-frame instant form, of classical rel- 
ativistic dynamics has been developed for those isolated systems of both N scalar particles 
with action-at-a-distance interactions and of charged scalar particles plus the electromag- 
netic field, which belong to timelike irreducible Poincare representations. While with only 
particles the kinematics is completely understood, in the case of particles plus fields more 
work will be needed to identify which point of the system has to be identified with the 
canonical 3-position of the center of mass on the spacelike hyperplane orthogonal to the 
total momentum, due to the presence of the interaction term in the constraints 7i p ~ of 
Eqs. (|146|) . This problem is non trivial and is connected with the action-reaction problem 
and with the classical basis of the electromagnetic interaction vertex of a scalar particle. In 
turn this point is associated with the pseudoclassical regularization of the classical electro- 
magnetic self-energies of charged particles. Again more work is needed to clarify the meaning 
of the pseudoclassical theory with Dirac observables where the effects of order ef are absent 
and only the e^-, i ^ j, effects are present: one has to compare the theory with Dirac 
observables to the standard methods of separation of the Coulomb and radiation effects in 
the Lienard-Wiechert potentials |55| , |56| , to revisit the Abraham-Lorentz-Dirac equation and 
the problem of classical bound states (Eqs.( |HS| ) define a regularized bound system with 
Coulomb interactions without emission of radiation and without fall on the centre). 

For free particles a quasi- Shanmugadhasan canonical basis has been found, which is 
suited to the description of the relativistic Cauchy problem for bound states, because the 
relative energies and times can be explicitely eliminated in a covariant way. It is under 
investigation [^6] how to introduce the second Poincare Casimir W 2 = —p 2 S of timelike 
Poincare representations in the basis and how to present in a clear form the gauge character 
of the classical zitterbewegung. The next step will be the introduction of the spin 1/2 degrees 
of freedom by means of Grassmann variables [0,0 both in the N- and 1-time theories and 
the comparison with the results coming from the Dirac equation. A further development will 
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be to replace the electromagnetic field with a non-Abelian Yang-Mills one and the search of 
the associated interparticle potential. 

When the electromagnetic field is considered, the 1-time theory shows clearly the exis- 
tence of a new formulation of the classical theory of free linear fields, i.e. the rest-frame field 
theory based on a center-of-mass and relative variables decomposition of field configurations. 
The formulation of rest-frame field theory, which is under active investigation, will be the 
necessary tool for attempting a quantization of the nonlocal field theories with Dirac's ob- 
servables, for utilizing the ultraviolet cutoff induced by the noncovariance of the canonical 
center-of-mass 3-position and for a fresh start to the problem of putting relativistic bound 
states among its Tomonaga-Schwinger-like asymptotic states. 

Furthermore, we have to understand the quantization of the 1-time covariant relativis- 
tic mechanics in the rest-frame instant form and to compare it with relativistic quantum 
mechanics [see the rich bibliography Refs. [|8T| , |S2| for the status of this theorywhich also 



originated from Ref. ||11|| 1, which is used in low energy (below the pion mass threshold) 
nuclear physics, where at the order 1/c 2 one can define a theory with a fixed number of 
nucleons. 

On the other hand, 1-time covariant relativistic statistical mechanics should be developed 
to a stage to be compared with the existing applications of relativistic statistical mechanics 



to relativistic (either astrophysical or nuclear) plasmas |79| 



In particular one should revisit the case of charged particles plus the electromagnetic field 
in the pseudo classical approach with Grassmann charges. One should obtain simplifications 
of the transport equations evaluated by taking into account the Abraham-Lorentz-Dirac 
equation with its dependence on higher accelerations More in general, it is still com- 



pletely open the problem of how to get a Hamiltonian formulation of the integro-differential 
equations of motion coming from the Tetrode-Fokker-Feynman- Wheeler actions, in which 
the radiation electromagnetic degrees of freedom have been eliminated: these actions admit 
reparametrization invariance of each particle worldline notwithstanding the interaction, are 
supposed to replace the nonexisting (due to the No-Interaction-Theorem) Lagrangians for 
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predictive mechanics [16| and are formally equivalent to an infinite system of differential 
eqations for accelerations of every order, whose study with constraint theory is just at the 
beginning p3| . 
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APPENDIX A: 



The rest frame form of the timelike fourvector p^ is V M = r\\fp 1 (\\ 0) = tf^spp 1 , P 2 = p 2 , 
where i] = signp°. 

o 

The standard Wigner boost transforming P M into p^ is 

o circ^ o 

\v ' o o — 

P p. V +p 2 

1 + 

z/ = e£(u(p)) = u^(p) = p^/r]\fp^ 

» = r e y(«(p)) = (-u r (p);y r - ^M ). (Al) 

o o 

The inverse of L M „(p,p) is L^ u (P,p), the standard boost to the rest frame, defined by 

= v(pi) = L"„(pi)|^- P - (A2) 

Therefore, we can define the following vierbeins [the e{f(u(p))'s are also called polarization 
vectors; the indices r, s will be used for A=l,2,3 and o for A=0] 

<£(u(p)) = L» A (p,°p) 

^(p)) = L\(p,p) = n AB r,^ B (u(p)) 



<£(«(p)) = V^ u (u(p)) = u^p) 

e r MP)) = -FVMp)) = (^(p);^-^ "^^ ) 
ef(u(p))=u A (p), (A3) 

which satisfy 

<f(«(p))<i(«(p)) = < 
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t A Mp)y B (u( P )) = ni 

rT = e» A (u(p)) V AB eUu(p)) = u»(p)u v {p) - £e{f(u(p)K(u(p)) 



r=l 



Vab = e^(u(p))^ 1/ e^(u(p)) 
The Wigner rotation corresponding to the Lorentz transformation A is 



R^(A,p) = [L(P,p)A- 1 L(Ap,P)} 



1 

i^(A,p) 



\ 



^•(A,p) = (A- 1 ) 1 , 



p P {A~ i y + vVp 1 

__ [{A -i )0j _ ((A-TW^A- 1 )^ 



(A4) 



(A5) 

p° + r]^ lK ^ >3 p p {k- 1 )P + r/v 7 ^ J ' 
The polarization vectors and the variables of Eq. (|13"D transform under the Poincare trans- 
formations (a, A) in the following way 



e?(u(Ap)) = (R-') r s A^e:(u(p)) 



= A^[x" + ^ rs i? r fc (A,p)^-i?\(A,p)] + = 
= AM y {y + - n Sr ; = [<(A° ,A ""- 1 ^ 



A° a p a + r\\fp 1 r s p° + ri\fp 1 
]} + a" 



(p u + VoVVp I )PrA° sl] /( 



r\\fp 1 {p° + r^i/p 7 

Tr' ar = ir as R s r (A,p). (A6) 

Therefore, x M is not a fourvector and p a , n a are Wigner spin-1 3- vectors; their infinitesimal 
transformation properties under Lorentz transformations generated by J^ u = L^ u + S^ v of 
Eq.(|4]), are 
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J af3 } = ri^x? - T]^x a + {x", S a/3 } 

s- } = 1 [jf r sl > + gfflg^ ] 

p° + Wp rjy'p^lp + 77 VP J 

= 

/n r TOM _ ^(/Pa-PX) 

S r T oi X $ is (p r <-p s <) 



P° + Tj^/p 1 

{tt^, J*} = (5 ts 5 jr - 6 ir 6 js )w s a . (A7) 

In contrast, the variables of Eqs. fl22|) transform under the Poincare transformations (a, A) 
in the following way 

T' = T + k fl {A- 1 af 
e = e 

z 'i = (A*. - £4- A ,) + e(A* - t-^A° )(A~ 1 a) At 
v J A ^ j; v M A° P P ^ A ; 



{T,p } = -v / l + P, {T,^} = -F 



{a* j oi } = -5 hi yi + p, {a* = 5 hj k i - 5 hi y 

1 + k 2 



{z h , J oi } = , {z h , J ij } = 5 hj z i - 8 hi z j . (A8) 

+ P 

Some further useful formulas are (e = 77VP 2 ) 

^-e B p {u{p)) = ^-L B p {p,p) = 
dp p p dpn 

= 2 ~foK - ^) - € 2 {p l +e / P B + ZVo)(P% + O + (p"tf + ^fop + 677°)] + 
(gB + £7? B)(p p + £7? °)[(p° + 2 e)p" + 6 2 ^] 
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1 d£(u(p)) 



pv 



- 2 VAC- 



I cAr r 

= --Wa(s* a - 

e UAK p° + e' 



S or p 



r p^ + 2erj 



e p B (u(p))S A « = 



e(p° + e) 



e(p° + e) 



(AlO) 



9 RMA,P) ^ 



dp 



( Pp (A-ry o + e y 



x 



x \-y + c(A- 1 )" )p j9 (A- 1 )^ - (p P (A _1 ) p + 6)(A~ 1 )^] - 
((A-^-l^A- 1 )^. p* 1 



r/ 



+ 



p° + e 
p° + e 



p p (A- 1 )" + e 



((A-To-l)x 



, „ 1 X + «T p" + e(A- 1 )% 



+ 



(p° + e) 2 e 



{-^ + e^)(A- 1 n + 



(AH) 
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APPENDIX B: 



In the N=2 case ]3§ Eqs.(|), (|), (0), 



' 5 VrM/ ' 



(^T|) become respectively 



P** = Pi + P2 



2^1 -P2); 



(Bl) 



;r 2 



Pi = ^ + Q 1 * 



i-i 

P2 



(B2) 



~v n , 1 a ( ( ^ de B {u{p)) 
= + - e v {u{p))rj AB — F — 5 p 



x r 



vVi^ip + vVp 2 ) 



dPv 



, PvP h 

p 2 



T, 



R 



p-R 

r}\[p I 

p-Q 



p r = R r - 

7T r = Q T - 



P 



r\yp 



p- R 



(Q°- 



p° + TJyjr' 

p-Q 



Tjy/p 2 p° + "q^pp 1 



(B3) 



1 



X- 



p-Q- \( m l ~ m l) = ee R - \( m i - m l) ~ 
1 



X+ = 2(0! + fc) = -(p 2 - M 2 + ){p 2 - M 2 _) + A X -(X- +mi-m 2 2 



P 
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i(e 2 - M 2 + )(e 2 - Ml) + 4 X -(X- + m\ - m\) « 



M± = ^mf - Qi ± ^/m 2 , -Q\ = \jm\ + vr 2 ± ^m 2 . + vr 2 , 

Q± = (»T - -jr)Q^ (B4) 











A" 


= i2" 







2p 2 [V p 2 ' 



1 2 

mi — m 



^FV, Q1 = Q1; (B5) 

z = x + g e ^ W))^ ^ p = 

= - 7^7 /as ^ + - sr^)) 

Vy/PVP+VvP) P 
% = T R 

mi — mi 

eR = e R x- = ee/? ~ 

2e 

P = P 

# = vr. (B6) 

Let us remark that Todorov's quasipotential-like Hamiltonian [44,4~Hj corresponds to the 
following rewriting of the constraint \+ (modulo \- ~ 0): 



X+ « 4[«ft + e p ^) 2 - m 2 ] = -4[tt 2 - 6 2 (e 2 )] « 



e 



m\m2 e 2 — m? — m 2 . 

m« = , e - 



; e A + m\ + ml-2{m\ + ml)e 2 -2m\ml ,, > 
b 2 (e 2 ) = i 2 ^| ^ ^ = t 2 p - m\- (B7) 



m p and e p are interpreted as the relativistic reduced mass and energy of a fictitious particle 
of relative motion. 
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APPENDIX C: 



Let {£(t)} be a one-parameter family of spacelike hypersurfaces foliating Minkowski 
spacetime M 4 . At fixed r, let z^(t, (?) be the coordinates of the points on E(r) in M 4 , {(?} 
a system of coordinates on S(r). If <r A = (<7 T = r;a = {<r r }) [the notation A = (r, f) with 
f = 1, 2, 3 will be used; note that A = r and A = f = 1, 2, 3 are Lorentz-scalar indices] and 
®A = d/da A , one can define the vierbeins 

4(r, a) = 9^(t, a), - d A z% = 0, (CI) 

so that the metric on S(r) is 

^(r, a) = -de* ||^(r, <?) II = (det\\z^(r, a) \ \f 

7(t,ct) = -det\\g f s{r,a)\\. (C2) 

If 7 rs (r, (?) is the inverse of the 3-metric g f s{T,a) [7™(r, <j)gus(r, (?) = the inverse 
# AB (r, <?) of c^(r, a) [# AC (r, a)^(r, a) = 5|] is given by 

7(t, <?) 



9 TT (r,a) = 



9(r, a) 



g»(r, (?) = Y s (t, (?) + [^7*7*1^ (C3) 



so that 1 = g rC {t, a)g^ T {r, (?) is equivalent to 



7(r, a) 

We have 



5- rr (r, a) - Y s (t, s ) 5Vr (r, 5) fiVs (r, <?) . ( C4) 



<M = d^ + g T ,Y s z^(r,a), (C5) 
V 7 



and 
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^ = z\{T,a)g AB {T,a)z v 6 {T,a) = 

= (n» + zh fS z u s )(r,a), (C6) 



where 



l 2 (r,a) = l, l tl (r,a)4(r,a) = 0, (C7) 

is the unit (future pointing) normal to S(r) at ^(r, cr). 
For the volume element in Minkowski spacetime we have 



<Pz = z?(t, a)drd% = dr«(r, a)Z M (r, a)) ^(r, a)d 3 a = 

lg{^f)drd?a. (C8) 



Let us remark that according to the geometrical approach of Ref . ]10| , by using Eq. (|C5 



in the form z£(r, a) = N{t, a)l^(r, a) + iV ? (r, 3)4 (t, a) [with N = and JV = 

being the standard lapse and shift functions, so that g TT = N 2 + gf S N r N s , g Tf = g r sN s , g TT = 
N~ 2 ,g Tf = -n f /N 2 ,g fS = 7" + we should write 

A T (r, a) = N(r, a)Mr, a) + N f (r, a)A f (r, a) (C9) 

and use Aj(r, a) as the genuine field configuration variable independent from the motion of 
the embedded hypersurface. 

Then, by denoting K n = K s ? = l^d f d- s z^ and rl = zffid^ = ~"f m (d?g sn + d s g ffl - 
dnOfs) the extrinsic curvature and the Christoffel symbols respectively of the spacelike hy- 
persurface (with metric g?$) embedded in the flat Minkowski spacetime, we get for the field 
strengths [";" and "|" are the total and ordinary covariant derivatives] 

= As 1 f — Af 1 § = dfA§ — Tf S An — (d§Af — T^A^) = dfA s — d§Af, 
F Tf = d T A f - df(NAi + N^Au) = NF lf + N S F S „ 

Fi? = A. f - A H = d f A x + K m ^ s A s -jj(dr- C-fjAf + \K r , r 'W, - A&N) = 
= ~[drAf - C$A f - fy(NAt)], (CIO) 



where C^A f = N S A?\ S + A s N s \ f = N s d s A f + A s d?N s is the Lie derivative along N. 

By using the second line of Eq.( p6|) to evaluate — \^grf v r\ pa F^ p F va 
-\N^y(2Y s Fi f Fi S + y^FfuFgv), the Lagrangian ( gig) becomes 

i N • • 

£ ( r^) = -Y,^-mm{r)ei{r)-9t{T)9 i {r)]- 
z i=i 

N i 

+ e^O^N^ a)Ai(r, a) + A f {r, a)(N f (r, 3) + tf(r))] - 



- \M r ' ^)[^7 rs («9 T ^ - CfiAf - dr(NAi))(d T A s - C^A S - d s (NA t )) + 

+ ^ T H 7 fifi J P«F„](r,?), (Cll) 

which is independent from both the extrinsic curvature and the Christoffel symbols [this is 
true only in the case of the spin 1 vector field but not for general tensors fields]. Now the 
configuration variables are z^(t, a), Ai(t, a), A?(t, a), and we perform the constraint analysis 
with the new associated momenta. However, since we have it 1 (t, a) = N(r,a)7i T (r,a) = 
and the other electromagnetic momenta still given by Eqs. (|115|) , in the case of (spin 1) 
vector fields both the formulation with Ai(t, a) and that with A t (t, a) give the same results 
as in Section VI, since both A t (t, a) and Ai(r, a) are gauge variables. 
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